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Theories with spontaneous local Lorentz and diffeomorphism violation contain massless Nambu- 
Goldstone modes, which arise as field excitations in the minimum of the symmetry-breaking poten- 
tial. If the shape of the potential also allows excitations above the minimum, then an alternative 
gravitational Higgs mechanism can occur in which massive modes involving the metric appear. The 
origin and basic properties of the massive modes are addressed in the general context involving an 
arbitrary tensor vacuum value. Special attention is given to the case of bumblebee models, which are 
gravitationally coupled vector theories with spontaneous local Lorentz and diffeomorphism violation. 
Mode expansions are presented in both local and spacetime frames, revealing the Nambu-Goldstone 
and massive modes via decomposition of the metric and bumblebee fields, and the associated sym- 
metry properties and gauge fixing are discussed. The class of bumblebee models with kinetic terms 
of the Maxwell form is used as a focus for more detailed study. The nature of the associated 
conservation laws and the interpretation as a candidate alternative to Einstein-Maxwell theory are 
investigated. Explicit examples involving smooth and Lagrange-multiplier potentials are studied to 
illustrate features of the massive modes, including their origin, nature, dispersion laws, and effects 
on gravitational interactions. In the weak static limit, the massive mode and Lagrange-multiplier 
fields are found to modify the Newton and Coulomb potentials. The nature and implications of 
these modifications are examined. 
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I. INTRODUCTION 

In relativistic quantum field theory, the nature of the 
field modes associated with the spontaneous breaking of 
an internal symmetry is now standard lore. When a 
global internal symmetry is spontaneously broken, one 
or more massless modes called Nambu-Goldstone (NG) 
modes appear If instead the symmetry is local, the 
Higgs mechanism can occur: the massless NG modes play 
the role of additional components of the gauge fields, 
which then propagate as massive modes [2]. In either 
case, the spontaneous symmetry breaking is typically 
driven by a potential term V in the Lagrange density. 
The vacuum field configuration lies in a minimum Vq of 
V. The massless NG modes can be understood as field 
excitations about the vacuum that preserve the value Vb, 
and they are associated with the broken generators of the 
symmetry. For many potentials, there are also additional 
excitations involving other values of V. These excita- 
tions, often called Higgs modes, correspond to additional 
massive modes that are distinct from any massive gauge 
fields. 

This standard picture changes when the spontaneous 
breaking involves a spacetime symmetry rather than an 
internal one. In this work, we focus on spontaneous 
breaking of Lorentz and diffeomorphism symmetries, for 
which the corresponding Higgs mechanisms exhibit some 
unique features (3|. Spontaneous Lorentz violation oc- 
curs when one or more Lorentz-nonsinglet field configu- 
rations acquire nonzero vacuum expectation values. The 
field configurations of interest can include fundamental 
vectors or tensors, derivatives of scalars and other fields, 
and Lorentz-nonsinglet composites. The nonzero vacuum 
values are manifest both on the spacetime manifold and 



in local frames Q ■ Their origin in spontaneous violation 
implies both local Lorentz violation and diffeomorphism 
violation, along with the existence of NG modes [f| . 

At the level of an underlying Planck-scale theory, nu- 
merous proposals exist that involve spontaneous Lorentz 
violation, including ones based on string theory [f| , non- 
commutative field theories Q, spacetime- varying fields 
Q, quantum gravity 0, random-dynamics models [l(|, 
multiverses [n|, brane- world scenarios 1121 . supersym- 
metry [l| , and massive gravity 0, 0, [l6[ . At experi- 
mentally accessible scales, the observable signals result- 
ing from Lorentz breaking can be described using ef- 
fective field theory [T3|. The general realistic effective 
field theory containing the Lagrange densities for both 
general relativity and the Standard Model along with 
all scalar terms involving operators for Lorentz violation 
is called the Standard-Model Extension (SME) 0, 
Searches for low-energy signals of Lorentz violation rep- 
resent a promising avenue of investi gati on involving the 
phenomenology of quantum gravity [19| . Numerous ex- 
perimental measurements of SME coefficients for Lorentz 
violation have already been performed 20], including 
ones with photon s I2li , electrons 
neutrons [25|, [2 a . [2 71 , mesons [28 
[3fj| . the Higgs [3l|7and gravity 

For spontaneous Lorentz and diffeomorphism break- 
ing, a general analysis of the nature of the NG modes 
and Higgs mechanisms is provided in Ref. [5j . One result 
is that the spontaneous breaking of local Lorentz sym- 
metry implies spontaneous breaking of diffeomorphism 
symmetry and vice versa. Since six local Lorentz trans- 
formations and four diffcomorphisms can be broken, up 
to ten NG modes can appear. To characterize these, it 
is natural to adopt the vierbein formalism [34| , in which 
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the roles of local Lorentz transformations and diffeomor- 
phisms are cleanly distinguished. It turns out that the 
vierbein itself naturally incorporates all ten modes. In an 
appropriate gauge, the six Lorentz NG modes appear in 
the antisymmetric components of the vierbein, while the 
four diffcomorphism NG modes appear along with the 
usual gravitational modes in the symmetric components. 

The dynamical behavior of the various NG modes is 
determined by the structure of the action 0]. In a La- 
grange density formed from tensor quantities and with 
diffeomorphism-covariant kinetic terms, the diffeomor- 
phism NG modes are nonpropagating. This feature, 
unique to diffcomorphism symmetry, can be viewed as 
arising because the diffeomorphism NG field excitations 
that preserve Vq include metric excitations, and the com- 
bined excitations cancel at propagation order in covari- 
ant derivatives and curvature. In contrast, the number 
of propagating Lorentz NG modes is strongly model de- 
pendent. For example, choosing the kinetic terms in the 
Lagrange density of the original field theory to eliminate 
possible ghost modes can also prevent the propagation of 
one or more Lorentz NG modes, leaving them instead as 
auxiliary fields. 

Several types of possible Higgs mechanisms can be dis- 
tinguished when spacetime symmetries are spontaneously 
broken. These have features distinct from the conven- 
tional Higgs mechanism of gauge field theories @ and 
Higgs mechanisms involving gravity without Lorentz vi- 
olation [35l | . The analysis of Higgs mechanisms can be 
performed either using the vierbein formalism, which per- 
mits tracking of Lorentz and diffeomorphism properties, 
or by working directly with fields on the spacetime man- 
ifold. The results of both approaches are equivalent. 

For local Lorentz symmetry, the role of the gauge fields 
in the vierbein formalism is played by the spin connec- 
tion. The Lorentz Higgs mechanism occurs when the 
Lorentz NG modes play the role of extra components 
of the spin connection [5( . Some components of the spin 
connection then acquire mass via the covariant deriva- 
tives in the kinetic part of the Lagrange density. Ex- 
plicit models displaying the Lorentz Higgs mechanism 
are known. For this mechanism to occur, the compo- 
nents of the spin connection must propagate as indepen- 
dent degrees of freedom, which requires a theory based on 
Riemann-Cartan geometry. If instead the theory is based 
on Riemann geometry, like General Relativity, the spin 
connection is fixed nonlinearly in terms of the vierbein 
and its derivatives. The presence of these derivatives en- 
sures that no mass terms emerge from the kinetic part of 
the Lagrange density, although the vierbein propagator 
is modified. 

In the context of diffcomorphism symmetry, the role of 
the gauge fields is played by the metric. For spontaneous 
diffcomorphism breaking in Riemann spacetime, a con- 
ventional diffeomorphism Higgs mechanism cannot gen- 
erate a mass for the graviton because the connection and 
hence the analogue of the usual mass term involve deriva- 
tives of the metric [3|. Also, since diffeomorphism NG 



modes are nonpropagating in a Lagrange density with 
covariant kinetic terms for reasons mentioned above, the 
propagating NG degrees of freedom required to generate 
massive fields via a conventional diffcomorphism Higgs 
mechanism are lacking. 

In a conventional gauge theory with a nonderivative 
potential V, the gauge fields are absent from V and so the 
potential cannot directly contribute to the gauge masses. 
However, in spontaneous Lorentz and diffeomorphism vi- 
olation, massive Higgs-type modes involving the vierbein 
can arise via the alternative Higgs mechanism, which in- 
volves the potential V 0. The key point is that V con- 
tains both tensor and metric fields, so field fluctuations 
about the vacuum value Vo can contain quadratic mass 
terms involving the metric. 

In this work, we study the nature and properties of the 
additional massive Higgs-type modes arising from this al- 
ternative Higgs mechanism. A general treatment is pro- 
vided for a variety of types of potentials V in gravita- 
tionally coupled theories with Riemann geometry. We 
investigate the effects of the massive modes on the phys- 
ical properties of gravity. In certain theories with spon- 
taneous Lorentz violation, the NG modes can play the 
role of photons 0] , and we also examine the effects of the 
massive modes on electrodynamics in this context. 

The next section of this work provides a general dis- 
cussion of the origin and basic properties of the massive 
Higgs-type modes. Section ILTTl analyzes the role of these 
modes in vector theories with spontaneous Lorentz vi- 
olation, known as bumblebee models. In Sec. IIVI the 
massive modes in a special class of bumblebee models 
are studied in more detail for several choices of poten- 
tial in the Lagrange density, and their effects on both 
the gravitational and electromagnetic interactions are ex- 
plored. Section fVl summarizes our results. Some details 
about transformation laws are provided in the Appendix. 
Throughout this work, the conventions and notations of 
Refs. [4|, [5| are used. 



II. MASSIVE MODES 

The characteristics of the massive Higgs-type modes 
that can arise from the alternative Higgs mechanism 
depend on several factors. Among them are the type 
of field configuration acquiring the vacuum expectation 
value and the form of the Lagrange density, including 
the choice of potential V inducing spontaneous break- 
ing of Lorentz and diffeomorphism symmetries. In this 
section, we outline some generic features associated with 
the alternative Higgs mechanism in a theory of a general 
tensor field T\^ u ... in a Riemann spacetime with metric 
9iw- We consider first consequences of the choice of po- 
tential V, then discuss properties of vacuum excitations, 
and finally offer comments on the massive modes arising 
from the alternative Higgs mechanism. 
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A. Potentials 

The potential V in the Lagrange density is taken to 
trigger a nonzero vacuum expectation value 

(T\fJ, V ...) = t\^y... (1) 

for the tensor field, thereby spontaneously breaking local 
Lorentz and diffeomorphism symmetries. In general, V 
varies with T\^ u ..., its covariant derivatives, the metric 
g^v, and possibly other fields. However, for simplicity 
we suppose here that V has no derivative couplings and 
involves only Tx^v- and g^. We also suppose that V 
is everywhere positive except at its degenerate minima, 
which have V = and are continuously connected via 
the broken Lorentz and diffeomorphism generators. The 
vacuum is chosen to be the particular minimum in which 
T\^ v ... attains its nonzero value |T]). 

Since the Lagrange density is an observer scalar, V 
must depend on fully contracted combinations of Tx^... 
and g^. Provided 7a m „... has finite rank, the number of 
independent scalar combinations is limited. For exam- 
ple, for a symmetric two-tensor field C^ v there are four 
independent possibilities, which can be given explicitly in 
terms of traces of powers of and g^ [36] . It is con- 
venient to denote generically these scalar combinations 
as X m , where to = 1, 2, . . . , M ranges over the number 
M of independent combinations. The functional form of 
the potential therefore takes the form 

V = V{X 1 ,X 2 ,...,X M ) (2) 

in terms of the scalars X m . 

The definition of the scalars X m can be chosen so that 
X m = in the vacuum for all m. For example, a choice 
involving a quadratic combination of the tensor with zero 
vacuum value is 

X = T x ^...g Xa g^g^ ■ ■ ■ T a01 ... ± t\ (3) 

with t 2 the norm 

f 2 = TV-(j A °>(/)(r)-^-, (4) 

where (g Xa ) is the vacuum value of the inverse met- 
ric. The =p sign is introduced for convenience so that 
t 2 can be chosen nonnegative. In principle, t 2 could vary 
with spacetime position, which would introduce explicit 
spacetime-symmetry breaking, but it suffices for present 
purposes to take t 2 as a real nonnegative constant. 

The M conditions X m = fix the vacuum value tx^v - 
If only a subset of N of these M conditions is generated 
in a given theory, then the value of t\p V ... is specified up 
to (M — N) coset transformations, and so the vacuum 
is degenerate under {M — N) additional gauge symme- 
tries. Note that these (M — N) freedoms are distinct 
from Lorentz and diffeomorphism transformations. 

It is useful to distinguish two classes of potentials V. 
The first consists of smooth functionals V of X m that are 
minimized by the conditions X m — for at least some 



to. These potentials therefore satisfy V — V r ' n — in the 
vacuum, where denotes the derivative with respect to 
X m , and they fix the vacuum value of Tx^v- to tx^v- 
modulo possible gauge freedoms. A simple example with 
the quantity X in Eq. (J3]) is 

V S (X) = \kX 2 , (5) 

where k is a constant. In this case, the vacuum value 
txfMu - is a solution of V = V = 0, where V denotes 
a derivative with respect to X. If the matrix V^ n of 
second derivatives has nonzero eigenvalues, the smooth 
functionals V can give rise to quadratic mass terms in the 
Lagrange density involving the tensor and metric fields. 
Potentials V of this type are therefore associated with the 
alternative Higgs mechanism, and they are the primary 
focus of our attention. 

A second class of potentials introduces Lagrange- 
multiplier fields X m for at least some to, to impose di- 
rectly the conditions X m = as constraint terms in 
the Lagrange density. We consider here both linear and 
quadratic functional forms for these constraints. An ex- 
plicit linear example is 

V L (X,X) = XX, (6) 

while a quadratic one is 

V Q (\,X) = \\X 2 . (7) 

In each example, A is a Lagrange-multiplier field and has 
equation of motion with solution X = 0, so the value 
txfj,u -. is a vacuum solution for the tensor. Potentials in 
the Lagrange-multiplier class arc unlikely to be physi- 
cal in detail because they enforce singular slicings in the 
phase space for the fields. However, when used with care 
they can be useful as limiting approximations to poten- 
tials in the smooth class, including those inducing the 
alternative Higgs mechanism [3]. Note that positivity of 
the potentials can constrain the range of the Lagrange 
multiplier field. For example, the off-shell value of A in 
Vl must have the same sign as X, while that of A in Vq 
must be non-negative. 

B. Excitations 

Field excitations about the vacuum solution {l} can be 
classified in five types: gauge modes, NG modes, massive 
modes, Lagrange-multiplier modes, or spectator modes. 
Gauge modes arise if the potential V fixes only (M — N) 
of the M conditions X m = 0, so that the vacuum is 
unspecified up to N conditions. These modes can be 
disregarded for most purposes here because they can be 
eliminated via gauge fixing without affecting the physics. 
The NG modes are generated by the virtual action of 
the broken Lorentz and diffeomorphism generators on 
the symmetry-breaking vacuum. They preserve the vac- 
uum condition V = 0. Massive modes are excitations for 
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which the symmetry breaking generates quadratic mass 
terms. Lagrange-multiplicr modes are excitations of the 
Lagrange multiplier field. Finally, spectator modes are 
the remaining modes in the theory. 

For smooth potentials, field excitations preserving the 
conditions X m — for all m have potential V — 0. They 
also satisfy = 0. The NG modes are of this type. Ex- 
citations with X m ^ having nonzero potential V =/= 
and 7^ are massive modes, with mass matrix deter- 
mined by the second derivatives V^ n . Since the smooth 
potentials depend on the tensor T\^ v ... and the metric 
g llVl the corresponding massive modes involve combina- 
tions of excitations of these fields. 

In contrast, for Lagrange- multiplier potentials the con- 
ditions X m = always hold on shell, which implies V = 
for all on-shell excitations. If also = 0, then the exci- 
tations remain in the potential minimum and include the 
NG modes. For linear functional forms of V, it follows 
that — A m , so V m is nonzero only when the Lagrange- 
multiplier field A m is excited. Any excitations of Txp,v- 
and <? M „ must have V = = 0. For quadratic functional 
forms of V, one finds = for all excitations, includ- 
ing the A field. We therefore can conclude that the com- 
binations of Tx^... and g M „ playing the role of massive 
modes for smooth potentials are constrained to zero for 
Lagrange-multiplicr potentials. Evidently, studies of the 
alternative Higgs mechanism must be approached with 
care when the Lagrange-multiplier approximation to a 
smooth potential is adopted. 

The tensor excitations about the vacuum can be ex- 
pressed by expanding Ta m „... as 



T 



(8) 



where the excitation tx^--- is defined as the difference 
txiiv- = 5Tx^... = T Aai „... - txnu- between the tensor 
and its vacuum value. We also expand the metric 



(9) 



in terms of the metric excitations h^v about the metric 
background value (g^)- For simplicity and definiteness, 
in much of what follows we take the background metric 
to be that of Minkowski spacetime, (<7 M „) — r/^. We also 
suppose tx^v - is constant in this background, so that 
d a txiw -- = in cartesian coordinates. 

Other choices can be made for the expansion of the 
tensor about its vacuum value. One alternative is to 
expand the contravariant version of the tensor as 



rj^XflU-- ■ £ 



XfJ.1 1 --- 



rpX^LV- ■ 



(10) 



In a Minkowski background, the vacuum values in the 
two expansions ([8]) and (fT0|) are related simply by 



t 



77 A< VV 7 ■ 



■t 



(ii) 



However, the relationship between the two tensor excita- 
tions at leading order involves also the metric excitation: 



rpX[lV- 



r \jJ,V- 



h Xa t j,u-_ hf ux t 



-h va t 



A// 



(12) 



In this expression, indices have been raised using the 
Minkowski metric. 

Any Lagrange multipliers A m in the theory can also be 
expanded about their vacuum values A m as 



A,; 



A m + A r 



(13) 



For linear Lagrange-multiplier potentials, the equations 
of motion for Tx^v -- and g^ v provide constraints on the 
vacuum values X m . In a Minkowski background and a 
potential yielding X m = for all m, the equations of 
motion for Tx^ v ... can be solved to give 



A m — 0. 



(14) 



For quadratic Lagrange-multiplicr potentials, the A m are 
absent from the equations of motion. Their vacuum val- 
ues are therefore physically irrelevant, and A m = can 
also be adopted in this case. For the remainder of this 
work we take A m = 0, and for notational simplicity we 
write A„, for both the full field A„, and the excitation A„, • 



C. Massive modes 

For a smooth potential V, the massive excitations typ- 
ically involve a mixture of tensor and gravitational fields. 
As an example, consider the expression for Vs in Eq. ([5]). 
This can be expanded in terms of the excitations tx^v--- 
and h^v. Retaining only terms up to quadratic order in 



a Minkowski background gives 



V S rj 2 K [^"'( V ..-i/ iA / )l ,..-i^t/ 



-i/l t 1 

2 W'A/i 



-...)] 2 , (15) 



where index contractions are performed with the 
Minkowski metric r]^ . Evidently, the massive excitations 
in this example involve linear combinations of ta m „... with 
contractions of and tx^.... 

The explicit expressions for the massive modes can be 
modified by local Lorentz and diffeomorphism gauge fix- 
ing. The action is symmetric under 10 local Lorentz 
and diffeomorphism symmetries, so there are 10 possi- 
ble gauge conditions. For an unbroken symmetry gener- 
ator, the effects of a gauge choice are conventional. For 
a broken symmetry generator, a gauge choice fixes the 
location in field space of the corresponding NG mode. 
For example, suitable gauge choices can place all the NG 
modes in the vierbein [5| • These gauge choices also affect 
the form of the massive modes. They can be used to iso- 
late some or all of the massive modes as components of 
either the gravitational or tensor fields by gauging other 
components to zero. Alternatively, the gauge freedom 
can be used to simplify the equations of motion, while 
the massive modes remain a mixture of the excitations 
TXft,,... and hpw 

The behavior of the massive modes depends on the 
form of the kinetic terms in the Lagrange density as well 
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as the form of the potential V. A gravitational theory 
with a dynamical tensor field T\^ v ... has kinetic terms for 
both g^ v and T\^ v ... and hence for both and t\^ v .... 
Since in the alternative Higgs mechanism the potential 
V acts only as a source of mass, the issue of whether the 
massive modes propagate dynamically depends on the 
structure of these kinetic terms. In particular, propa- 
gating massive modes can be expected only if the theory 
without the potential V contains the corresponding prop- 
agating massless modes. 

It is desirable that any propagating modes be unitary 
and ghost free. To avoid unitarity issues with higher- 
derivative propagation, the kinetic term for the metric 
excitation h^ v can be taken to emerge as usual from the 
Einstein-Hilbert action. In the absence of the potential 
V, this is also a ghost-free choice. For the tensor field, 
higher-derivative propagation can be avoided by writ- 
ing the general kinetic term Ck as a weighted sum of 
all scalar densities formed from quadratic expressions in 
T\fi U ... that involve two covariant derivatives. For exam- 
ple, one such scalar density is 

C K ~ eT x ^...D a D a T x ^-, (16) 

where e = \J—g is the vierbein determinant. In the 
absence of the potential V, the ghost-free requirement 
places strong constraints on allowed forms of £k and typ- 
ically involves gauge symmetry for the tensor. For a vec- 
tor field, for example, the Maxwell action is the unique 
ghost-free combination in the Minkowski spacetimc limit. 
In any case, since the massive modes are combinations of 
h^v and tx m „..., it follows that ghost-free propagation is 
possible only if the kinetic terms for these combinations 
are ghost free. Note that the potential V may explicitly 
break the tensor gauge symmetry, so requiring ghost-free 
kinetic terms is by itself insufficient to ensure ghost-free 
massive modes. 

Under the assumption of Lorentz invariance, the Fierz- 
Pauli action involving quadratic terms for is the 
unique ghost-free choice for a free massive spin-2 field 
[37| . However, when coupled to a covariantly conserved 
energy-momentum tensor, the small-mass limit of a mas- 
sive spin-2 field includes modes that modify the gravita- 
tional bending of light in disagreement with observation 
[381 ]. This presents an obstacle to constructing a viable 
theory of massive gravity. One avenue of investigation 
that might permit evading this obstacle is to relax the 
assumption of Lorentz invariance. Spontaneous Lorentz 
violation in closed string theory has been proposed [l4| 
as a mechanism that might lead to graviton mass terms. 
Models involving infrared modifications of gravity have 
also been proposed that involve spontaneous Lorentz vi- 
olation with ghosts [l5| and that have explicit Lorentz 
violations [1 61 ] . Explicit Lorentz violation is incompatible 
with Ricmann and Riemann-Cartan geometries but may 
be compatible with Finsler or other geometries [3, [3£fl 
or may be viewed as an approximation to spontaneous 
violation. In the present context, the possibility exists 
that the massive modes from the alternative Higgs mech- 



anism could evade the Veltman-van Dam-Zakharov con- 
straint via their origin in spontaneous Lorentz violation 
and their nature as mixtures of gravitational and ten- 
sor modes. For example, the expansion (|15[) includes 
quadratic terms involving contractions of with the 
vacuum value t\^ v , so a model with suitable vacuum val- 
ues and incorporating also a ghost-free propagator for 
the corresponding modes would describe a type of prop- 
agating massive gravity without conventional Fierz-Pauli 
terms. 

In the alternative Higgs mechanism, the existence of 
propagating massive modes involving the metric can be 
expected to affect gravitational physics. Effects can arise 
directly from the modified graviton propagator and also 
from the massive modes acting as sources for gravita- 
tional interactions. The latter can be understood by 
considering the energy-momentum tensor T^ v of the full 
theory, which can be found by variation of the action 
with respect to the metric g^. The contribution Ty V to 
arising from the potential V is 

T^ = -g^V + 2V^^, (17) 

where a sum on m is understood. For the vacuum solu- 
tion X m — 0, which satisfies V = V 7 ' n = 0, Ty V remains 
zero and the gravitational sourcing is unaffected. The 
same is true for excitations for which V and both van- 
ish, such as the NG modes. However, the massive modes 
arising from a smooth potential have nonzero V and Vj n 
and can therefore act as additional sources for gravity. 
These contributions can lead to a variety of gravitational 
effects including, for example, modifications of the New- 
ton gravitational potential in the weak-field limit, which 
could have relevance for dark matter, or cosmological fea- 
tures such as dark energy 

Note that nonpropagating modes can also have sim- 
ilar significant physical effects on gravitational proper- 
ties. This possibility holds for any excitations appear- 
ing in V, whether they are physical, ghost, or Lagrange- 
multiplier modes. For example, any Lagrange-multiplier 
fields are auxiliary by construction and so the A m ex- 
citations are nonpropagating. Nonetheless, for linear 
Lagrange-multiplier potentials these excitations can con- 
tribute to Ty V even though V = because is nonva- 
nishing. However, a theory with a quadratic Lagrange- 
multiplier potential has Ty U = and therefore leaves the 
gravitational sourcing unaffected. 

III. BUMBLEBEE MODELS 

In this section, we focus attention on the special class of 
theories in which the role of the tensor T\^... is played by 
a vector that acquires a nonzero vacuum expectation 
value b^. These theories, called bumblebee models, are 
among the simplest examples of field theories with spon- 
taneous Lorentz and diffcomorphism breaking. In what 
follows, bumblebee models are defined, their properties 
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under local Lorentz and diffeomorphism transformations 
are presented, their mode content is analyzed, and issues 
involving gauge fixing and alternative mode expansions 
are considered. 



A. Basics 

The action Sb for a single bumblebee field coupled 
to gravity and matter can be written as 

Sb = J d A x Cb 

= J d A x (C g +C gB +C K +Cv+C M )- (18) 

In Riemann spacetime, the pure gravitational piece C g of 
the Lagrange density is usually taken to be the Einstein- 
Hilbert term supplemented by the cosmological constant 
A. The gravity-bumblebee couplings are described by 
C g B , while £k contains the bumblebee kinetic terms and 
any self-interaction terms. The component Cy consists 
of the potential V(B^), including terms triggering the 
spontaneous Lorentz violation. Finally, Cm involves the 
bumblebee coupling to the matter or other sectors in the 
model. 

The forms of C 9 b and Ck are complicated in the 
general case. However, if attention is limited to terms 
quadratic in B^ involving no more than two derivatives, 
then only five possibilities exist. The Lagrange density 
Cb can then be written as 

C B = -—e{R-2\) + <j 1 eB»B v R tiv + <y 2 eB>*B li R 
-IneB^B^ + \r- 1 eD ll B v B^B v 
^\T Z eD^D v B v -eV + £ M , (19) 

where G is the Newton gravitational constant and the 
field-strength tensor B^ v is defined as 

= d^B v - d u B„, (20) 

The five real parameters ci, a 2 , 72, 13 in Eq. (fTO)) 
are not all independent. Up to surface terms, which leave 
unaffected the equations of motion from the action, the 
condition 

J d 4 x {eB»B v R^ u - \eB» v B, w 

+eD fl B l/ D^B' J - eD^DvB") = (21) 

is identically satisfied, so only four of the correspond- 
ing five terms in Cb are independent. No term of the 
form eB^B^R appears in the condition (f2"Tj) , so the pa- 
rameter <72 remains unaffected while the four parameters 
cii T ij T 2, 73 become linked. For example, the identity 
(l2"Tj) implies that the action for the Lagrange density (flU)) 
with five nonzero parameters <7i, a 2 , 71, r%, T3 is equiv- 
alent to an action of the same form but with only four 



nonzero parameters <r[ — o\ + 73, a' 2 — c%, t[ = t\ — 2t^, 
t' 2 = T2 — 2T3, while T3 = 0. Moreover, other factors may 
also constrain some of the five parameters. For example, 
certain models of the form Cb yield equations of motion 
that imply additional relationships among the parame- 
ters. Also, certain physical limits such as the restric- 
tion to Minkowski spacetime can limit the applicability 
of the condition (f2"Tj) . Some cases with specific param- 
eter values may be of particular interest for reasons of 
physics, geometry, or simplicity, such as the models with 
Ti = 1, T2 = T3 = considered below, or the model with 
02 = — ci/2 for which the bumblebee-curvature coupling 
involves the Einstein tensor, Cb D <j\eB^ B v G ^ . Since 
the most convenient choice of parameters depends on the 
specific model and on the physics being addressed, it is 
useful to maintain the five-parameter form (|19[) for gen- 
erality. 

Following the discussion in Sec. Ill Al we suppose the 
potential V in Eq. (fT!))) has no derivative couplings and is 
formed from scalar combinations of the bumblebee field 
B^ and the metric . Only one independent scalar 
exists. It can be taken as the bumblebee version of X in 
Eq. ©: 

X = B^ V B V ± b 2 , (22) 

where b 2 is a real nonnegative constant. The potential 
V(X) itself can be smooth in X, like the form ([5]), or 
it can involve Lagrange multipliers like the form ([6]) or 
([?])■ In any case, the vacuum is determined by the single 
condition 

X = Bf,g^B u ± b 2 = 0. (23) 

In the vacuum, the potential vanishes, V{X) = 0, and 
the fields B^, acquire vacuum values 

-» (B^) = b„, 5m „ -> (g^). (24) 

The nonzero value 6 M , which obeys b fJi (g tl ' y )b L , = =F& 2 , 
spontaneously breaks both Lorentz and diffeomorphism 
symmetry. Note that the choice of the potential V can 
also have implications for the parameters in Eq. (|19p . For 
example, if the potential takes the quadratic Lagrange- 
multiplier form |7[). then <ii can be taken as zero be- 
cause a nonzero value merely acts to rescale G and A. 
However, a nonzero value of <ii can have nontrivial con- 
sequences for models with other potentials, such as the 
smooth quadratic form ([5]). 

For generic choices of parameters, the Lagrange den- 
sity (|19p is unitary because no more than two derivatives 
appear. However, as discussed in Sec. Ill CI the indefinite 
metric and generic absence of gauge invariance typically 
imply the presence of ghosts and corresponding negative- 
energy problems, which can tightly constrain the viabil- 
ity of various models ■ If the gravitational couplings 
and the potential V are disregarded, there is a unique 
set of parameters ensuring the absence of ghosts: n = 1, 
T2 = T3 = 0. With this choice, the kinetic term for 
the bumblebee becomes the Maxwell action, in which 
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the usual U(l) gauge invariance excludes ghosts. When 
the gravitational terms and couplings are included, this 
gauge invariance is maintained for the parameter choice 
T\ = 1, Ta = T3 = o\ = (T2 = 0. The further inclusion of 
the potential V breaks the U(l) gauge symmetry, but the 
form of the kinetic term ensures a remnant constraint on 
the equations of motion arising from the identity 

D f j l D v B> J,v = 0. (25) 

The action for these bumblebee models, introduced in 
Ref. [3|], is therefore of particular interest. The corre- 
sponding Lagrange density 

£kS = J^Q< R - 2K )-\ eB ^ B ^- eV + C m, (26) 

is investigated in more detail in Sec. IIVI The reader is 
warned that some confusion about the relationship be- 
tween these models and ones with nonzero T2 and T3 ex- 
ists in the literature. In particular, results for the mod- 
els (|2T)1) can differ from those obtained in models with 
nonzero T2, T3 via straightforward adoption of the limit 
T 2,T3 — * 0, due to the emergence of the remnant con- 
straint (ggj) . 

The discussion above defines bumblebee models on the 
spacetime manifold, without introducing a local Lorentz 
basis. In this approach, the Lorentz NG modes are hid- 
den within the bumblebee field B^. Adopting instead a 
vierbein formulation reveals explicitly the local Lorentz 
properties of the models, and it also provides a natural 
way to incorporate spinor fields in the matter Lagrange 
density Cm ■ The vierbein e^ a converts tensors expressed 
in a local basis to ones on the spacetime manifold. The 
spacetime metric <? M „ is related to the Minkowski metric 
rjab hi the local frame as 

g»v = e^ej'rjab, (27) 

while the bumblebee spacetime vector B^ is related to 
the local bumblebee vector B a as 

= e^ a B a . (28) 

A complete treatment in the vierbein formalism involves 
also introducing the spin connection w^" 6 , which appears 
in covariant derivatives acting on local quantities. In 
Riemann-Cartan geometry, where the spacetime has both 
curvature and torsion, the spin connection represents de- 
grees of freedom independent of the vierbein. However, 
experimental constraints on torsion are tight (4lj . In 
this work, we restrict attention to Riemann geometry, 
for which the torsion vanishes and the spin connection 
is fixed in terms of the vierbein. It therefore suffices 
for our purposes to consider the vierbein degrees of free- 
dom and their role relative to the gravitational and NG 
modes. Bumblebee models in the more general context 
of Riemann-Cartan spacetime with nonzero torsion are 
investigated in Refs. 0, [I[ . 



There is a substantial literature concerning theories 
of vacuum- valued vectors coupled to gravity. The five- 
parameter Lagrange density (fT9)) excluding the poten- 
tial V and the cosmological constant was investigated by 
Will and Nordtvedt in the context of vector-tensor mod- 
els of gravity [52]. Kostelecky and Samuel (KS) [H in- 
troduced the potential V triggering spontaneous Lorentz 
violation and studied both the smooth quadratic poten- 
tial ([5]) and the linear Lagrange-multiplier case ([6]) for 
the class of models given by Eq. (|26[) . The presence of 
the potential introduces a variety of qualitatively new 
effects, including the necessary breaking of U(l) gauge 
invariance the existence of massless NG modes and 
massive modes [Bj], and potentially observable novel ef- 
fects for post-newtonian physics [33j and for the matter 
sector [43J. The potential V also leads to a candidate 
alternative description of the photon [5j and the gravi- 
ton 0, H(| • In Minkowski spacetime, more general po- 
tentials V of hypergeometric form are known to satisfy 
the one-loop exact renormalization group [Hj]. Models 
of the form (|19|) with <j\ — 02 — 0, a unit timelike 6 M , 
a linear Lagrange-multiplier potential Vl, and an addi- 
tional fourth-order interaction for £? M have been studied 
in some detail as possible unconventional theories of grav- 
ity 0, |4(|. Other works involving bumblebee models 
include Refs. 0, HH M, M, [EH . 

An aspect of bumblebee models of particular interest 
is the appearance of massless propagating vector modes. 
This feature suggests the prospect of an alternative to the 
usual description of photons via U(l) gauge theory. The 
central idea is to identify the photon modes with the NG 
modes arising from spontaneous Lorentz violation. Early 
discussions along these lines centered on reinterpretating 
the photon or electron in the context of special relativ- 
ity without physical Lorentz violation (52l. |53||. More re- 
cently, the Lorentz NG modes in certain bumblebee mod- 
els with physical Lorentz violation have been shown to 
obey the Einstein-Maxwell equations in Riemann space- 
time in axial gauge Q . These models are further consid- 
ered below, with the discussion initiated in Sec. IIV A 31 
One motivation of the present work is to investigate the 
role of massive modes and Lagrange-multiplier fields in 
this context. In particular, it is of interest to investi- 
gate possible modifications to electrodynamics and grav- 
ity, with an eye to novel phenomcnological applications 
of spontaneous Lorentz and diffcomorphism breaking. 



B. Transformations 

In considering theories with violations of Lorentz and 
diffcomorphism symmetry, it is important to distinguish 
between observer and particle transformations [181 ]. Un- 
der either an observer general coordinate transformation 
or an observer local Lorentz transformation, geometric 
quantities such as scalars, vectors, tensors, and their 
derivatives remain unchanged, but the coordinate basis 
defining their components transforms. In contrast, parti- 
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cle transformations can change geometric quantities, in- 
dependently of any coordinate system and basis. 

In theories without spacetime-symmetry breaking, the 
component forms of the transformation laws for particle 
and observer transformations are inversely related. For 
example, under infinitesimal particle diffcomorphisms de- 
scribed by four infinitesimal dispacements the com- 
ponents of the bumblebee field transform according to 
the Lie derivative as 

-> - [d M ?)B v - nB„, 
B» -> B* + {d v e)B v - Z v d v B», (29) 

while the metric transforms as 

g» v -> g» v - (df,C)gau - (dv£ a )g»a - Cd a g^, 

g*» _ g ^ + {d a e)g av + (9 Q f)r - rv. (30) 

Under infinitesimal observer general coordinate transfor- 
mations, which are the observer equivalent of diffcomor- 
phisms, the formulae for the transformations of the bum- 
blebee and metric components take the same mathemat- 
ical form up to a possible sign change in the arbitrary pa- 
rameter even though these transformations are only 
the result of a change of coordinates. Similarly, under 
infinitesimal particle local Lorentz transformations with 
six parameters e a b = ~£ba related to the local Lorentz 
group element by A a b 5 a b + e a b , the local components 
of the bumblebee field transform as 

B a -> B a + e a b B b , (31) 

with the formula for B a following from this by raising in- 
dices with rj ah . Under observer local Lorentz transforma- 
tions, the corresponding transformation formulae again 
take the same mathematical form. However, any form of 
physical Lorentz and diffeomorphism breaking destroys 
the details of these equivalences. 

A fundamental premise in classical physics is that the 
properties of a physical system are independent of the 
presence of a noninteracting observer. The observer 
is free to select a coordinate system to describe the 
physics of the system, but the physics cannot depend 
on this choice. In particular, this must remain true even 
when Lorentz and diffeomorphism symmetry are broken, 
whether explicitly or spontaneously. Viable candidate 
theories with spacetime-symmetry breaking must there- 
fore be invariant under observer transformations, which 
are merely changes of coordinate system. For example, 
the SME is formulated as a general effective field the- 
ory that is invariant under observer general coordinate 
transformations and under observer local Lorentz trans- 
formations H 111 El- 

In contrast, a theory with physical Lorentz and dif- 
feomorphism breaking cannot by definition remain fully 
invariant under the corresponding particle transforma- 
tions. For example, if the breaking is explicit, the action 
of the theory changes under particle transformations. If 
instead the breaking is spontaneous, the action remains 



invariant and the equations of motion transform covari- 
antly, as usual. However, the vacuum solution to the 
equations of motion necessarily contains quantities with 
spacetime indices that are unchanged by particle trans- 
formations. These vacuum values and the excitations 
around them can lead to physical effects revealing the 
symmetry breaking. 

In bumblebee models, the relevant spacetime vacuum 
values are those of the bumblebee and metric fields, de- 
noted (B,j,) and (g^v). These are unaffected by particle 
diffeomorphisms, and their nonzero components thereby 
reveal the broken particle diffcomorphisms. In contrast, 
(Bfj) and (g^J) both transform as usual under observer 
general coordinate transformations, which therefore are 
unbroken. Analogous results hold for the local-frame vac- 
uum value (B a ) of the local bumblebee field. The com- 
ponents (B a ) remain unaffected by particle local Lorentz 
transformations, with the invariance of the nonzero com- 
ponents resulting from the breaking of local Lorentz gen- 
erators, while the components (B a ) transform under ob- 
server local Lorentz transformations in the usual way. 
Similarly, the vacuum value (e„ a ) of the vierbein is un- 
changed by both particle diffeomorphisms and particle 
Lorentz transformations, but it transforms as a vector 
under the corresponding observer transformations. 

By virtue of the relation Q28p between the spacetime 
and local bumblebee fields, it follows that spontaneous lo- 
cal Lorentz violation is necessarily accompanied by spon- 
taneous diffeomorphism violation and vice versa [5[ . The 
point is that the vacuum value of the vierbein is nonzero, 
so the existence of a nonzero (B a ) spontaneously break- 
ing local Lorentz symmetry also implies the existence of 
a nonzero (B^) spontaneously breaking diffeomorphism 
symmetry, and vice versa. Note, however, that this re- 
sult fails for explicit violation, where the analogue of the 
relation (|28|) is absent. In general, explicit local Lorentz 
violation occurs when a nonzero quantity t a b c ... is exter- 
nally prescribed in the local frame, but the corresponding 
spacetime quantity tx^v... = e^ a e„ e„ c • • • t a &c... is defined 
using the full vierbein and hence remains invariant un- 
der diffeomorphisms. Similarly, explicit diffeomorphism 
violation occurs when a nonzero quantity t\^ v _. is exter- 
nally prescribed on the spacetime manifold, but the cor- 
responding local quantity defined via the inverse vierbein 
remains invariant under local Lorentz transformations. 

C. Mode expansions 

To study the content and behavior of the modes, the 
fields can be expanded as infinitesimal excitations about 
their vacuum values. At the level of the Lagrange den- 
sity, it suffices for most purposes to keep only terms to 
second order in the field excitations, which linearizes the 
equations of motion. 

Assuming for simplicity a Minkowski background 
(gnu) — i)y,vi we expand the metric and its inverse as 

V = V + V- ff"" w t?"" - ft"". (32) 
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For the local bumblebee vector, we write 

B a = b a +(3 a , (33) 

where the vacuum value (B a ) is denoted b a and assumed 
constant, and the infinitesimal excitations are denoted 
f3 a . The vacuum condition (|23p implies 

baV ab h = Tb 2 . (34) 

In a Minkowski background, the vierbein vacuum value 
can be chosen to be (e^ a ) = 8^ in cartesian coordinates. 
The expansions of the vierbein and its inverse are 

e'a » ^-|^ + /«. (35) 

where /i Ma and x^a are, respectively, the symmetric and 
antisymmetric components of the vierbein. The covari- 
ant and contravariant components of the bumblebee field 
follow from Eq. (gSJ) as 

B, « S^(b a + f3 a ) + (lh^ + X , a )b a , 

fl" * ^ a ( 6 « + j g«) + (_ i^ + ^jft-. ( 36 ) 

Since the local and spacetime background metrics are 
both Minkowski and since the excitations are infinitesi- 
mal, the distinction between the Latin local indices and 
the Greek spacetime indices can be dropped. We adopt 
Greek indices for most purposes that follow, raising and 
lowering indices on purely first-order quantities with rf' v 
and rifiu- The vierbein expansions (|35|) can then be 
rewritten as 



total of 20. On the right-hand side, the symmetric met- 
ric component has 10, the antisymmetric component 
Xnv has 6, the transverse bumblebee excitation has 
3, and the longitudinal excitation (3 has 1, producing the 
required matching total of 20. Of these 20 degrees of free- 
dom, 6 are metric modes, 4 are bumblebee modes, while 
6 are associated with local Lorentz transformations and 
4 with diffeomorphisms. The explicit transformations of 
all the field components under particle diffeomorphisms 
and local Lorentz transformations can be deduced from 
the full-field expressions and from the invariance of the 
vacuum values. A list of formulae for both particle and 
observer transformations is provided in the Appendix. 

Among all the excitations, only the longitudinal com- 
ponent (3 of the bumblebee field is invariant under both 
particle diffeomorphisms and local Lorentz transforma- 
tions. It is therefore a physical degree of freedom in any 
gauge. Moreover, using Eqs. ([23| and (j40|) reveals that 
exciting the [3 mode alone produces a nonzero value of 
X, given at first order by 

X = B^ V B V ± b 2 » 2(6%)/?. (41) 

As a result, the excitation (3 is associated with a nonmin- 
imal value of the potential, and it therefore cannot be an 
NG mode. In fact, for the case of a smooth quadratic 
potential, (3 is a massive mode. In contrast, in a theory 
with a Lagrange- multiplier potential where the constraint 
X = is enforced as an equation of motion, the massive 
mode (3 identically vanishes. 

D. Gauge fixing and NG modes 



« rf a ^ a w if - \h» v + x^- (37) 

The vacuum value for the bumblebee vector becomes 

(Bp) = (e/>6 - V'6a = V (38) 

It is convenient to decompose the local vector excita- 
tions /3fj, = 5^ a (3 a into transverse and longitudinal pieces 
with respect to b^. Excluding for simplicity the case of 
lightlight bfj,, we write 

/3 M = /3*+/36„, /?> M = 0, (39) 

where b^ = b^/Vlr* is a vector along the direction of 
bfj, obeying 6^& M = =pl. Using this decomposition, the 
bumblebee mode expansions (|36p become 

B„ w V + (iV+X„/+^ + 4 

B 11 sa b" + + x llv )K+ (3^ + Pb". (40) 

It is instructive to count degrees of freedom in the ex- 
pressions (|3"T)) and (|4"U|) . On the left-hand side, the vier- 
bein has 16 components and the bumblebee field 4, for a 



Since the spacetime-symmetry breaking in bumble- 
bee models is spontaneous, the four diffeomorphisms 
parametrized by £^ and the six local Lorentz transforma- 
tions parametrized by e^u leave invariant the bumblebee 
action (|19p and transform covariantly the equations of 
motion. Fixing the gauge freedom therefore requires 10 
gauge conditions. 

For the diffeomorphism freedom, a choice common in 
the literature is the harmonic gauge 

= 0, (42) 

where = h^ u — ^rj^h and h = —h = —h^^. In 
the harmonic gauge, the Einstein tensor becomes as 
— \\^h^ v at linear order. An alternative choice for the 
diffeomorphism degrees of freedom is the axial gravita- 
tional gauge 

h^b" = 0. (43) 

Both these gauge choices represent four conditions. To 
fix the local Lorentz freedom, it is common to eliminate 
all six antisymmetric vierbein components by imposing 
the six conditions 

Xp» = 0- (44) 
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Other choices are possible here too. Consider, for ex- 
ample, the decomposition of Xuv m terms of projections 
along by, 

Xfxv — Xfiu ~F Xffiv i Xu^fii 
X%V = x'y = 0, (45) 

which is the analogue of Eq. (|3"9")) for f3y. The compo- 
nents Xyv an d xji = XinM 1 each contain three degrees 
of freedom. Inspection of the transformation laws shows 
that an alternative to fixing the local Lorentz gauge via 
Eq. (|44|) is the set of six conditions 

X% =01=0- (46) 

Note, however, that the combination (xyub v + fly) ap- 
pearing in By in Eq. (|40[) is invariant and therefore can- 
not be gauged to zero. Evidently, the associated local 
Lorentz degrees of freedom must remain somewhere in 
the theory. 

The bumblebee vacuum value (By) = by breaks one 
of the four diffeomorphism symmetries. The broken gen- 
erator is associated with the projected component £ v b u 
of the parameter Analogously, the vacuum value 
(B a ) = b a breaks three of the six local Lorentz sym- 
metries. The broken generators are associated with the 
components e a bb b Ri e.y U b u of ty V projected along b v . In 
each case, the unbroken generators are associated with 
the orthogonal complements to the projections. 

Since the vacuum breaks one particle diffeomorphism 
and three local Lorentz transformations, four NG modes 
appear. A useful general procedure to identify these 
modes is first to make virtual particle transformations 
using the broken generators acting on the appropriate 
vacuum values for the fields in V, and then to promote 
the corresponding parameters ty V and £ M to field excita- 
tions [f|: 

The properties of S M and £ M „ under various particle and 
observer transformations are given in the Appendix. The 
projections ^ v b v and ty V b v associated with the broken 
generators determine the NG modes, which are therefore 
E v b v and S^V. 

We can follow this procedure to elucidate the relation- 
ship between the NG modes and the component fields 
in the decomposition (|40|) of By. Consider first the dif- 
feomorphism NG mode, which is generated by a virtual 
transformation acting on the spacetime vacuum value 
{By) and involving the broken diffeomorphism genera- 
tor. The relevant transformation is given in Eq. (|All) . 
and it yields 

(B^ -^by- (d^E v )b v . (48) 

Comparison of this result to the form of Eq. (1401) reveals 
that the four vierbein combinations (^h^ + Xtiv)b v con- 
tain the diffeomorphism NG mode. This agrees with the 



result obtained by combining virtual diffeomorphisms on 
the vacuum values of the component fields in Eq. (|4"0"1) : 

(Pi) -> 0, (J3) - 0. (49) 

Note also that the combinations {khn V +Xnv)W maintain 
the potential minimum V = 0, as is expected for an NG 
mode. 

In contrast, the three Lorentz NG modes £y V b v are 
generated by virtual transformations acting on the local- 
frame vacuum value (B a ) and involving the broken local 
Lorentz generators: 

(B a ) ^b a +£ ab b b , (50) 

Comparison with Eqs. (|33p and (|39|) shows that the 
transverse field /?' contains the three Lorentz NG modes 
£tu/b u . Note that there are exactly three degrees of free- 
dom in (3* , all of which maintain the potential minimum 
V = 0. This result can be used to connect the com- 
ponent fields in the decomposition of By with the three 
Lorentz NG modes £y V b v . The vacuum value (By) itself 
is invariant under local Lorentz transformations, so an- 
other component field must also contain the Lorentz NG 
modes. Performing virtual local Lorentz transformations 
on the vacuum values of the component fields gives 

(VO -> o, ( X „vb u ) -> -£ MV r, 

((3D -> E^V, (13) -» 0, (51) 

which shows that the three combinations (xyv^ v ) also 
contain the three Lorentz NG modes. 

Combining the results (|4"9"]) and (|5ip reveals the follow- 
ing NG mode content for the component field combina- 
tions in the decomposition (|40| of By. 

{-\hy V + x^vW = -(dyB v )b v - E^b", 
S3 1 — F b v 

t'y — °fJ,vU ! 

D 0. (52) 

We see that the four combinations (^h^ + X/j.u)b v con- 
tain a mixture of the diffeomorphism NG mode and the 
three Lorentz NG modes, while the three combinations 
/3* contain only the three Lorentz NG modes. 

By fixing the 10 gauge freedoms in the theory, the 
above results can be used to determine the physical con- 
tent of the field B^. Suppose first we adopt the 10 condi- 
tions (|4"3"1) and (|4"4")) . Then, the bumblebee field becomes 

By » by + Pl + (3b^ 

= by + £yX + /3by, (53) 

where the result ([52]) has been used. Note that the fixed 
gauge means that fields with different transformation 
properties can appear on the left- and right-hand sides. 
In this gauge, the four components of By are decomposed 
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into three Lorentz NG modes associated with ft. and the 
one massive mode (3. The diffeomorphism NG mode is 
absent. It obeys 

{d„E v )V = -{d v E^)b v (54) 

and hence b fl (d fl E u )b L ' — 0, and it is locked to the Lorentz 
NG modes via the equation 

(flpS„)6" = -£„„&". (55) 

For the diffeomorphism NG mode, this gauge is evidently 
analogous to the unitary gauge in a nonabelian gauge 
theory. 

An alternative gauge choice could be to impose the 10 
conditions (|43|) and (|46|) instead. This gives 

« b^ + xl + Ph 

= b^ + S^ + pb^, (56) 

where in this gauge the Lorentz NG modes E^b" are 
identified with x^ instead of /3* . One way to understand 
this identification is to perform a local Lorentz transfor- 
mation with parameter e^„ = —E^u on the first result 
in Eq. ([53")) . This changes the value of ft from E^J? 
to zero, while simultaneously converting from zero to 
Xjt = E^yW . In this gauge, the explicit decomposition 
(l56|) of Bp in terms of NG modes is the same as that 
in Eq. (|53p . but the three Lorentz NG modes are now 
associated with the components of the vierbein. The 
diffeomorphism mode remains absent. It still obeys the 
condition ([54)) and is locked to the Lorentz NG modes by 
Eq. (E9 . 

Partial gauge conditions can also be imposed. For 
example, suppose the choice (|4*51) is made for the local 
Lorentz gauge, while the diffeomorphism gauge remains 
unfixed. Then, the four degrees of freedom in the com- 
bination i^h^ v + Xtiv)b v consist of the three Lorentz NG 
modes and the diffeomorphism NG mode, and all the NG 
modes are contained in the vierbein [5j. The bumblebee 
field can be written as 

w ^ - {d^E v )b u + E^b v + $b M . (57) 

Note that only one projection of E v appears, even though 
four diffeomorphism choices remain. The corresponding 
expression for B^ includes additional metric contribu- 
tions and is given by 

B" m W + {d u E»)b v + E^b u + ftb". (58) 

This equation contains contributions from all four fields 
E^. However, if the diffeomorphism excitations are re- 
stricted only to the one for the broken generator, for 
which S M obeys b^E u = b v E^, then also reduces 
to an expression depending only on the diffeomorphism 
NG mode E v b v . Related results are obtained in Ref. 
Q , which investigates the fate of the NG modes using a 
decomposition of the vierbein into transverse and longi- 
tudinal components along b^. This decomposition leads 



to the same relations as Eqs. ([57)) and ((58)) when the 
condition b^E v — b v E^ is applied. 

Even without gauge fixing, the fields E^ cancel at lin- 
ear order in both and B^. As a result, propa- 
gating diffeomorphism NG modes cannot appear. This 
is a special case of a more general result. By virtue 
of their origin as virtual particle transformations, the 
diffeomorphism NG modes appear as certain compo- 
nents of representation-irreducible fields with nonzero 
vacuum values. However, diffeomorphism invariance en- 
sures these modes enter in the metric and bumblebee 
fields in combinations that cancel in a diffeomorphism- 
invariant action, including the general bumblebee action 
(flU)) . In contrast, the Lorentz NG modes do play a role 
in the bumblebee action. They can be identified with a 
massless vector field = E^V in the fixed axial gauge 
sa 0. The fully gauge-fixed expression (|53| for B^ 
then becomes 

B^^b^ + A^+ftb^ (59) 

where the transverse components have the form of photon 
fields in the axial gauge, and the longitudinal mode is the 
massive mode (3. 

E. Alternative expansions 

The analysis in the previous subsections is based on 
the local- frame expansion ([3"5")) of the bumblebee field B a . 
However, other mode expansions are possible, including 
ones in which the vierbein makes no explicit appearance 
and the local Lorentz transformations are no longer man- 
ifest. This subsection offers a few comments on two al- 
ternative expansions used in some of the literature. 

The first alternative mode expansion is specified using 
the covariant spacctime components of the bumblebee 
field 

B fl = b fl + E fi , (60) 

where E^ represents the excitations of the spacetime 
bumblebee field around the vacuum value b^. It follows 
that the contravariant components are 

B^^ + ^-ri),,. (61) 

These fields are linked to the vierbein and the local-frame 
fluctuations via 

e M = $h lu , + x»v)V + 0l+0b li - (62) 

The fields E^ are scalars under particle local Lorentz 
transformations, but transform under particle diffeomor- 
phisms as 

- - (d^ a )b a , 

E» = ^ V E V -> £» -if v {d v t, a )b a . (63) 
Note that this usage of £ M differs from that in Ref. [j| . 
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The second alternative expansion (33l.l45l.l47ll4i.l49t 

starts instead with the contravariant bumblebee compo- 
nents 

B t * = b t *+£'», (64) 

where a prime is used to distinguish the excitations £' M 
from the previous case. The corresponding covariant 
components are then 

B lt f*b lt + 8' ll + h ltu V'. (65) 

The field redefinitions connecting these to the vierbein 
and previous case are 

= { r \hf iv + 7t u ')b v +ffi + f3b" 

= ifE^-hTK. (66) 

The fields are scalars under particle local Lorentz 
transformations but transform under particle diffcomor- 
phisms as 

£'V _> £ '» + (d a ^)b a , 
£'„ = v^E" 1 ^ £'„ + (8 a ^)b a . (67) 

Note that the fields £ M and £' have different transforma- 
tion properties and contain different mixes of the bum- 
blebee and metric excitations. However, £^ and £'^ take 
the same gauge- fixed form in the gauge (|4*3"|) . 

The two alternative expansions are useful because the 
excitations £ M and £ ' are invariant under infinitesimal lo- 
cal particle Lorentz transformations. In these variables, 
the local Lorentz symmetry is hidden and only the dif- 
feomorphism symmetry is manifest. The Lorentz NG 
modes lie in the 6^-tranverse components of £ M or £'^. 
The massive mode joins the diffeomorphism mode in 
the longitudinal component of £ M . It can be identified 
from Eq. (|62[) as the combination 

& = ^iziV^ = Tb , { ^ _ i h ^ n (68) 

b a b a 

Notice that (3 appears here as a diffeomorphism-invariant 
combination of field components that transform nontriv- 
ially. 

The choice of diffeomorphism gauge has interesting 
consequences for £^ and £'. First, note that £^ can- 
not be gauged to zero, since only one degree of freedom 
£ a b a appears in its transformation law. However, it is 
possible to gauge £'^ to zero. The corresponding gauge 
condition for £ M is £^ = h^b" . In either of these gauges, 
the massive mode becomes 

P = h^ir, (69) 

thereby becoming part of the metric. The bumblebee 
field strength reduces to 



and so is also given by the metric. Evidently, in these 
gauges the theory is defined entirely in terms of the met- 
ric excitations h^ v . Requiring £^ to vanish or dropping 
the term h^b" in Eq. (f6"5)) therefore improperly sets to 
zero the bumblebee field strength B^, which alters the 
equations of motion governing the Lorentz NG and mas- 
sive modes [HI]. It follows from Eq. d70j) that the ef- 
fective Lagrange density can contain additional kinetic 
terms for h^, beyond those arising from the gravita- 
tional action, that originate from the bumblebee kinetic 
terms £k- These additional terms propagate the Lorentz 
NG modes, disguised as massless modes contained in the 
metric. 



IV. EXAMPLES 

In this section, some simple examples are developed to 
illustrate and enrich the general results obtained in the 
discussions above. We choose to work within the class 
of KS bumblebee models with Lagrange density given by 
Eq. (|26p . which avoids a priori propagating ghost fields 
and the complications of nonminimal gravitational cou- 
plings. For definiteness, we set A = and choose the 
matter-bumblebee coupling to be 

C M = -eB^J M , (71) 

where J M is understood to be a current formed from mat- 
ter fields in the theory. In principle, this current could 
be formed from dynamical fields in the theory or could 
be prescribed externally. For the latter case, diffeomor- 
phism invariance is explicitly broken unless the current 
satisfies a suitable differential constraint. 

The section begins by presenting some general results 
for this class of models, including ones related to the 
equations of motion, conservation laws, and the con- 
nection to Einstein-Maxwell electrodynamics. We then 
turn to a more detailed discussion of three specific cases 
with different explicit potentials V(X), where the bum- 
blebee field combination X is defined as in Eq. (|22[) . 
The three cases involve the smooth quadratic potential 
Vs(X) of Eq. ((SJ), the linear Lagrange-multiplier poten- 
tial Vl(A, X) of Eq. (J6]), and the quadratic Lagrange- 
multiplier potential Vq(X,X) of Eq. ([7]). The massive 
modes are studied for each potential, and their effects on 
gravity and electromagnetism are explored. 



A. General considerations 

1. Equations of motion and conservation laws 

Varying the Lagrange density (|26|) with respect to the 
metric yields the gravitational equations of motion 



Bfj, u = (d^h va — d v h,j,„)b a , 



(70) 



= 8ttGT^. 



(72) 
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In this expression, the total energy-momentum tensor 
T^ v can be written as a sum of two terms, 

T"" = T&" +Tg. (73) 

The first is the energy-momentum tensor TjT' arising 
from the matter sector. The second is the energy- 
momentum tensor T^ v arising from the bumblebee ki- 
netic and potential terms, 

where and Ty V are given by 

= B» a B\ - y^B aP B a ? (75) 

and 

T$ v = -Vg» y + 2V I B^B U . (76) 

Varying instead with respect to the bumblebee field 
generates the remaining equations of motion, 

D V B^ = J*. (77) 

Like the total energy-momentum tensor, the total current 
J p can be written as the sum of two terms 

J» = J& + J%- (78) 

The partial current is associated with the matter sec- 
tor and acts as an external source for the bumblebee field. 
The partial current arises from the bumblebee self- 
interaction, and it is given explicitly as 

J£ = -2V'B». (79) 

The contracted Bianchi identities for Gy V lead to con- 
servation of the total energy-momentum tensor, 

DyT^ = D^Ttf + Tg") = 0. (80) 

The antisymmetry of the bumblebee field strength By V 
implies the remnant constraint 125)1 an d hence a second 
conservation law, 

= D„{J^ + j£) = 0. (81) 

Note that this second conservation law is a special fea- 
ture of KS bumblebee models. It is a direct consequence 
of choosing the a priori ghost-free action (|26p. Note also 
that this conservation law holds even though the poten- 
tial term V excludes any local U(l) gauge symmetry in 
these models. 



2. Bumblebee currents 

The bumblebee current Jg defined in Eq. (|79p van- 
ishes when V' = 0. This situation holds for the vac- 
uum solution and for NG modes, and it then follows from 



Eq. I)81p that the matter current is covariantly con- 
served. However, V' =/= in the presence of a nonzero 
massive mode or a nonzero linear Lagrange multiplier, 
whereupon the bumblebee current Jg acts as an addi- 
tional source for the bumblebee field equation (|77|) . A 
nonzero massive mode or nonzero linear Lagrange mul- 
tiplier contributes to the bumblebee component Tg" of 
the energy-momentum tensor too, and it therefore also 
acts as an additional source for the gravitational field 
equations . 

The contributions to the energy-momentum tensor Tg" 
stemming from V ^ can be positive or negative. This 
is a generic feature, as can be seen from the expression 
for Ty V in Eq. (17) and the general bumblebee Lagrange 
density (|T9|) . The full energy-momentum tensor T^ v is 
conserved, but the prospect of negative values for 
implies stability issues for the models. Under suitable 
circumstances, for example, unbounded negative values 
of Tg can act as unlimited sources of energy for the 
matter sector. Stability is also a potential issue in the 
absence of matter. For example, for the case where is 
disregarded and the potential involves a linear Lagrange 
multiplier producing A ^ on shell, at least one set of 
initial values is known that yields a negative-energy so- 
lution 

Whether instabilities occur in practice depends on the 
form of the Lagrange density and on the initial condi- 
tions. The situation is comparatively favorable for KS 
bumblebee models because the extra conservation law 
(|8ip can play a role. In principle, negative-energy contri- 
butions in the initial state can be eliminated by choosing 
initial conditions such that V = 0, while the conser- 
vation law can prevent the eventual development of a 
destabilizing mode. One way to see this is to expand the 
conservation law Eq. (|8ip at leading order in a Minkowski 
background. For nonzero (Bo), we obtain 

doV « ^(« - 2V'd li B» - 2B>d 3 V'). (82) 

Suppose the matter current is independently conserved, 
DyJ^ ss dyJ^ = 0. The initial condition V = then 
implies 8qV' = initially. Taking further derivatives 
shows that all time derivatives of V vanish initially and 
so V remains zero for all time, indicating stability is 
maintained. 

It is interesting to note that the matter-current conser- 
vation law D^jjf — emerges naturally in the limit of 
vanishing NG and massive modes, where the bumblebee 
field reduces to its vacuum value B^ — > (By) = by. This 
agrees with a general conjecture made in Ref. [4j. The 
line of reasoning is as follows. In the limit, the bumble- 
bee theory with spontaneous Lorentz violation takes the 
form of a theory with explicit Lorentz violation with cou- 
plings only to by. However, theories with explicit Lorentz 
violation contain an incompatibility between the Bianchi 
identities and the energy-momentum conservation law. 
The incompatibility is avoided for spontaneous Lorentz 
violation by the vanishing of a particular variation in the 
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action, which reduces in the limit to the requirement of 
current conservation D^jjf = 0. Nonetheless, if is 
excited away from this limit, current conservation in the 
matter sector may fail to hold. 

In practice, the potential instability may be irrelevant 
for physics. For example, if a bumblebee model is viewed 
as an effective field theory emerging from an underlying 
stable quantum theory of gravity, any apparent insta- 
bilities may reflect incomplete knowledge of constraints 
on the massive modes or of countering effects that come 
into play at energy scales above those of the effective the- 
ory. Under these circumstances, in KS bumblebee models 
it may suffice for practical purposes simply to postulate 
that the matter and bumblebee currents do not mix and 
hence obey separate conservation laws, 

D^J M = 0, D^ B = 0. (83) 

For instance, one can impose that only matter terms Cm 
with a global U(l) symmetry are allowed, as is done for 
the case of Minkowski spacetime in Ref. [50| . Another 
option might be to disregard Jjf altogether and allow 
couplings between the matter and bumblebee fields only 
through gravity. In what follows, we adopt Eqs. (|83"|) 
as needed to investigate modifications of gravity arising 
from massive modes and to study the conditions under 
which Einstein-Maxwell solutions emerge in bumblebee 
models. 



3. Bumblebee electrodynamics 

An interesting aspect of the KS bumblebee models is 
their prospective interpretation as alternative theories of 
electrodynamics with physical signatures of Lorentz vio- 
lation |5j. In this approach, basic electrodynamics prop- 
erties such as the masslessness of photons are viewed as 
consequences of spontaneous local Lorentz and diffeomor- 
phism breaking rather than of exact local U(l) symmetry. 
A key point is that KS bumblebee models have no depen- 
dence on second time derivatives of Bq. They therefore 
feature an additional primary constraint relative to the 
bumblebee models (|19[) with more general kinetic terms. 
The primary constraint generates a secondary constraint 
in the form of a modified Gauss law, which permits a 
physical interpretation of the theory in parallel to elec- 
trodynamics. This modified Gauss law is unavailable to 
other bumblebee models. 

In a suitable limit, the equations of motion (|72p and 
(177)) for the KS bumblebee models reduce to the Einstein- 
Maxwell equations in a fixed gauge for the metric and 
photon fields. To demonstrate this, we work in asymp- 
totically Minkowski spacetime and choose the gauge in 
which the bumblebee field is given by Eq. ([55]). The limit 
of interest is /5 — s- 0, corresponding to zero massive mode. 
For the case with the Lagrange-multiplier potential ©, 
where /3 is constrained to zero, an equivalent result also 
follows in the limit A — > 0. Since Ty V acquires nonzero 



contributions only from the massive mode or linear La- 
grange multiplier, it follows that as/3— » or A — > the 
bumblebee energy-momentum tensor Tg" reduces to that 
of Einstein-Maxwell electrodynamics, 

Tg" - IT = = F^FoT - b^F aP F^. (84) 

In this equation, the bumblebee field strength is reinter- 
preted as the electromagnetic field strength 

B lu , = F^, = d ll A v -d v A^ (85) 

via the mode expansion (|59[) . and the gravitational field 
equations (|72|) reduce to their Einstein-Maxwell counter- 
parts. Similarly, the bumblebee current Jg in Eq. (|79|) 
vanishes as (3 — * or A — » because V' — > 0. This 
leaves only the conventional matter current J M , which 
by virtue of Eq. (f8Tj) obeys covariant conservation. The 
bumblebee field equations f77| therefore also reduce to 
their Einstein-Maxwell counterparts. 

Even when (3 vanishes, the interpretation of the model 
as bumblebee electrodynamics can in principle be distin- 
guished from conventional Einstein-Maxwell electrody- 
namics through nontrivial SME-type couplings involving 
the vacuum value b^ in the matter equations of motion. 
Depending on the types of couplings appearing in the 
matter-sector Lagrange density Cm, a variety of effects 
could arise. At the level of the minimal SME, for ex- 
ample, axial vector couplings involving can be sought 
in numerous experiments, includin g on es with electrons 
[H[23[, protons and neutrons [25Ll26j|. and muons [29j ] . 
Other types of SME coefficients can also be generated 
from the vacuum value b^. For example, the nonzero 
symmetric traceless SME coefficients might emerge 
via the symmetric traceless product & M 6„ — 6 a 6 Q! 7y Att/ /4, 
and experimental searches for these include ones with 
electrons [12, H3] , protons, and neutrons [H, Hfij ■ 

If the model conditions permit the field (3 to be 
nonzero, further deviations from Einstein-Maxwell elec- 
trodynamics can arise. For example, in the weak-field 
regime where (3b ^ is small compared to b lll the bumblebee 
current Jg is linear in (3 and acts as an additional source 
in the equation of motion (|77|) for the field strength B^ v . 
Effects of this type are investigated in the next few sub- 
sections. For large fields, nonlinear effects in the bum- 
blebee current and energy-momentum tensor may also 
play a role. However, the field value of j3 is typically 
suppressed when the mass of (3 is large, so if the mass is 
set by the Planck scale then the linear approximation is 
likely to suffice for most realistic applications. 

In the limit of zero f3, the physical fields are the field 
strengths i? M „ and the gravitational field. Excitations of 
the bumblebee field B^ in the classical theory are then 
unobservable, in parallel with the potential A^ in clas- 
sical electromagnetism. This may have consequences for 
leading-order effects in the weak-field limit. For exam- 
ple, the weak static limit in general bumblebee models 
(fT9|) requires B^ itself to be time independent. However, 
at leading order in the KS bumblebee model, it can suf- 
fice to require that only the field strength B^ v be time 
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independent while £? M itself has time dependence, in anal- 
ogy with Maxwell electrodynamics in Minkowski space- 
time. For example, the Coulomb electric field E = — 8qA 
for a static point charge q emerges in temporal gauge 
A = from a time-dependent — (0,tdj& q ), where 
$ q is the Coulomb potential. Likewise, in certain leading- 
order weak static solutions for which j3 acts as a source 
of charge, the bumblebee field may naturally exhibit 
potential lines converging at the source of charge, similar 
to the physical configurations with singular derivatives of 
A^ that occur in classical electrodynamics. 

Another interesting issue is the role of quantum ef- 
fects. Since bumblebee electrodynamics involves gravity, 
it faces the same quantization challenges as other theo- 
ries of gravity and electrodynamics, including Einstein- 
Maxwell theory. Some work on the renormalization 
structure at one loop has been performed 0, [55[ , and in 
the limit of zero massive mode and in Minkowski space- 
time the usual properties of quantum electrodynamics 
are expected to hold. Addressing this issue in detail is 
an interesting open problem that lies beyond the scope 
of this work. 



B. Smooth quadratic potential 

This subsection considers the specific KS bumblebee 
model ([25]) having smooth quadratic potential §5§ with 
the definition (JUJ): 

V = V S (X) = ± K (B^B„±b 2 ) 2 . (86) 

For definitcness, we adopt the mode expansion B^ = 
bn + £n of Eq. (|60p in a Minkowski background and as- 
sume weak fields and £^ so that the gravitational 
and bumblebee equations of motion (j?2")) and (|77|) can be 
linearized. 

One goal is to investigate deviations from Einstein- 
Maxwell theory arising from the presence of a weak 
nonzero massive mode (3. We therefore focus on domi- 
nant corrections to the linearized Einstein-Maxwell equa- 
tions arising from (3. In this approximation, the bumble- 
bee component of the energy-momentum tensor becomes 

T b ~ T em + Mb a b a W l b»P, (87) 

where is the zero-/3 limit of given in Eq. ([84]) and 
the other term arises from Ty . The bumblebee current 
reduces to 

J£ n -4K(b a b a )b"/3. (88) 

The constraint obtained by assuming current conserva- 
tion in the matter sector becomes 

6^/3 w 0. (89) 

These expressions reveal that at leading order the mas- 
sive mode (3 acts as a source for gravitation and electro- 
dynamics, subject to the constraint (|89|) . 



Note that the linearization procedure can alter the 
time behavior and dynamics of the fields in the presence 
of a nonzero massive mode (3. Suppose, for example, that 
6^ is timclike and we adopt the global observer frame in 
which fop = (6,0,0,0). Nonlinearities in the current 
associated with nonzero (3 generate time dependence for 
most solutions because the spatial current Jb does work 
on the field strength B^. However, the linearization 
(1551 implies Jg is nonzero only along the direction of 
& M , reducing to a pure charge density Jg ~ [ps, 0,0,0) 
with pb ~ —4Kb 2 P, and the current-conservation law 
(|5"9"| then requires pb and hence (3 to be static. Nonethe- 
less, the linearization procedure captures the dominant 
effects of nonzero (3. 

1. Propagating modes 

To study free propagation of the gravitational and 
bumblebee fields in the absence of charge and matter, 
set = = in the linearized equations of motion. 
The gravitational equations (|72[) then become 

Dh^u + v^d a d h a /3 - d a d l Ji av - d a dJh aVL 

w -MTrGK{b a b a )b^b v (3. (90) 

Contributions from are second-order in £ M and can 
be neglected in this context. The bumblebee equations 
([77ll reduce to 

□ £„ - d„d u £ v » Mb a b a )b^p. (91) 

Note that the massive mode (3 is defined by Eq. (j6"5|) in 
terms of £^ and h^, so it is not an independent field in 
these equations. Note also that the linearized current- 
conservation law (|89p follows by taking a derivative of 
Eq. (jS). 

To investigate the behavior of the massive mode j3, it is 
convenient to combine the above two equations to obtain 

uP-^{b a b a ){i+AKGbPbp)[3 « -hr-&dx& t (e li -h MV r). 

b'<b 1 

(92) 

At first glance, this equation might suggest that (3 can 
be a propagating massive field. However, j3 depends on 
the fields appearing on the right-hand side, so the modes 
are still coupled and care is required in determining the 
dispersion properties, including the mass value. 

A suitable choice of diffcomorphism gauge clarifies 
matters. The modes remain entangled in the harmonic 
gauge (|42|) and also in the barred axial gravitational 
gauge h^W = 0. However, a sufficient decoupling can be 
achieved by adopting the axial gravitational gauge (14"3")1 
and by decomposing £ M into pieces parallel and perpen- 
dicular to b^. In this gauge, the longitudinal component 
£ of £fj, reduces to £ = (3, so that 

£^ = £l + £b^ = A fi +^b fl . (93) 
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In this equation, we write £^ = for the transverse 
components of £^ to make easier the task of tracking /in- 
dependent deviations from conventional electrodynamics. 
These components satisfy the axial condition = 0. 

We emphasize that the axial condition is not a gauge 
choice. It is a consequence of projecting along 6 M and is 
independent of the gravitational gauge fixing. It does, 
however, have the same mathematical form as an axial- 
gauge condition in electrodynamics, even though the 
bumblebee models have no U(l) symmetry. 

With these choices, the usual Einstein equations for 
h^ v in axial gravitational gauge are recovered from the 
result (|90j) when (3 vanishes Q. For nonzero f), only one 
linearly independent combination of the Einstein equa- 
tions is changed. It can be written 

Uh + 2d^d u h llv « ~MTiGn{b a b a ){b p bp)(3. (94) 

Since the usual propagation-transverse components of 
/i M „ remain unaffected, the physical graviton modes prop- 
agate normally. The bumblebee equations (|9Tj) become 

□ A" — d^d^A v + y^-b fl b a d a d v A v w (95) 

and 

UP - Mb a b a )f3 « -^b a d a d»A^ (96) 
bPbp 

When the massive mode (3 vanishes, these equations re- 
duce to those of electrodynamics in axial gauge A 
nonzero value of (3 modifies the equations, but the usual 
propagation-transverse components of A^ and hence the 
physical photon modes remain unaffected. 

Since f3 and A^ are independent fields, the dispersion 
properties of the massive mode (3 can be identified from 
Eq. (|96|) . Any solutions of the theory describing freely 
propagating modes are required to obey the boundary 
conditions that the spacetime be asymptotically flat and 
that the bumblebee fields vanish at infinity. At linear 
order, these solutions are formed from harmonic plane 
waves with energy-momentum vectors p^ — {E,p) obey- 
ing suitable dispersion relations. The modes satisfying 
the asymptotic boundary conditions are then constructed 
as Fourier superpositions of the plane- wave solutions. For 
the massive mode f3, the constraint equation (|89|) imposes 
the additional requirement 

&% « 0. (97) 

Any freely propagating massive mode is therefore con- 
strained to have an energy-momentum vector orthogonal 
to the vacuum value 6 M . For harmonic- mode solutions 
to the equations of motion, this requirement implies that 
the term on the right-hand side of Eq. (|96p vanishes in 
Fourier space. The resulting dispersion law for the mas- 
sive mode /3 is then 

p% « -Anb a b a , (98) 



which involves the squared-mass parameter 

Mf 3 = Anb a b a . (99) 

The sign of M§ depends on whether b^ is timelike or 
spacelike. The scale of Mj depends on both k and b. 
Inspection of Eqs. (f9"4")) or (fM]) shows that the limits 
IM^I 2 — > oo or Hi — ► oo with 6 M fixed are equivalent to 
taking the limit of vanishing massive mode, (3 — ► 0. The 
discussion in Sec. IIV A 31 then implies that the limit of 
large |M^| 2 approximates Einstein-Maxwell theory. 

For the case of a timelike vacuum value b^, the pa- 
rameter Mg — —4Kb 2 has the wrong sign for a physi- 
cal mass. Adopting the global observer frame in which 
b^ = (b, 0, 0, 0), we see that the constraint (|9T|) forces the 
energy E of this mode to vanish, leaving the condition 
that the magnitude of the momentum must remain fixed 
at \p\ = \Mp\. A time- independent mode with fixed spa- 
tial wavelength cannot be Fourier superposed to form a 
physical wave packet satisfying the asymptotic bound- 
ary conditions. It follows that no physical propagating 
massive mode appears when b^ is timelike. 

If instead the vacuum value b^ is spacelike, then Mj = 

4k& 2 is positive and is a candidate for a physical mass. 
Choose the global observer coordinate system such that 
the spacelike 6 M takes the form 6 M = (0, 0, 0, b). The con- 
straint (|97p then requires the vanishing of the momentum 
along the spatial direction 63, and the dispersion law ([98]) 
becomes 

E 2 -(p 2 x +p 2 y )ttM 2 g . (100) 

This suggests that (3 could propagate as a harmonic plane 
wave along spatial directions transverse to 6 M . How- 
ever, these harmonic plane waves are constant in z and 
hence fail to satisfy the asymptotic boundary conditions 
as z — > ±00 unless their amplitude is zero. Since all 
harmonics propagate under the same constraint (|97[) . no 
Fourier superposition can be constructed to obey the 
boundary conditions. It follows that no physical prop- 
agating massive mode exists for the case of a spacelike 
vector either. 

The above results can also be confirmed by obtaining 
the eigenvalues of the Fourier-transformed equations of 
motion, which form a 14x14 matrix determining the na- 
ture of the modes. Consider, for example, the case of 
purely timelike b^ . In the axial gravitational gauge, four 
of the modes are zero, corresponding to the four fixed 
degrees of freedom. Another four eigenvalues describe 
modes propagating via a massless dispersion law, cor- 
responding to the two photon and two graviton modes. 
One eigenvalue has zero energy and sets |p| = \Mp\, cor- 
responding to the massive mode (3. The remaining 5 
modes are auxiliary, and all have zero energy. Related 
results can be obtained for the case of spacelike 6 M . 
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2. Weak static limit 

Although the massive mode (3 is nonpropagating, it 
plays the role of an auxiliary field and can thereby pro- 
duce various effects. For example, its presence can affect 
the weak static limits of the gravitational and bumblebee 
interactions. In particular, a nonzero (3 can modify the 
forms of both the Newton and Coulomb potentials. 

To demonstrate this, we consider the weak static limit 
of the field equations with an external matter sector con- 
taining massive point charges. The matter Lagrange den- 
sity can be taken as 



dr dr 



1/2 



S 4 (x - x n (r)) 



n 

+ J2ln / dr B^x)^S\x-x n (r)), (101) 

where m n and q n are the masses and charges of the point 
particles, respectively. The energy-momentum tensor for 
the matter sector is 

t m = \ m « / dr 



dx^ L dx 
dr dr 



n r4 



5\x-x n (r)). (102) 



Its combination with Tg given in Eq. (|74| forms the 
total energy-momentum tensor, which is conserved. The 

{Pqi J) is 



matter-sector current = 



7^ 



dr^S 4 (x 
dr 



(r)). 



(103) 



0. 



It is conserved, D^J^ 

The standard route to adopting the weak static limit 
is to linearize the equations of motion in the fields and 
then discard all time derivatives. This meets no diffi- 
culty under suitable circumstances, such as the vacuum 
value being purely spacelike, and we follow this route 
in the present subsection. However, in certain cases ad- 
ditional care is required because the 6^-transverse com- 
ponents of the bumblebee field obey the axial condition 
6 p £* t = 0, which can imply time dependence even for 
static fields. For example, as discussed in Sec. IIV A 31 
the case of purely timclikc and zero (3 is equivalent to 
electrodynamics in temporal gauge, for which the static 
Coulomb potential involves a linear time dependence in 
A. The corresponding weak static limit with nonzero (3 
can therefore be expected to exhibit this feature. We 
return to this issue in the next subsection. 

To proceed, it is useful to choose a diffeomorphism 
gauge that simplifies the Einstein tensor and allows di- 
rect comparison between the geodesic equation and the 
Newton force law. One possible choice is the harmonic 
gauge (j4"2"|) , in which the static gravitational potential Q g 
is related to the metric by $ s = — ^h no = — tj/ioo 
The linearized gravitational equations become 



where f3 is the massive mode given in terms of £ M and 
h^ u in Eq. (JSHJ), and is the zero-/? form of T£ v given 
in Eq. ([El]). 

We can now extract the relevant equations in the weak 
static limit. The equation for the gravitational potential 
$ s is the 00 component of the trace-reversed form of Eq. 

cm 



v 2 $„ 



87rG[T£? + T°° M + 2 K (b a b a )(b 2 + b 2 )(3], (105) 



where the trace-reversed energy-momentum tensors are 
defined by T^ v = — rj^T a a and are evaluated in 
the static limit. The weak static bumblebee equations 
become 



V 2 £ 

V 2 £- V(V-f) 



Pq + 4 K (b a b a )b (3, 
J + 4k(6"S q )6/3, 



(106) 



huh 



EM 



4k(& q 6 q )^6 1/ /3], (104) 



where p q and J are static charge and current sources. 

The above equations show that the massive mode (3 
can be viewed in the weak static limit as a simultane- 
ous source of energy density, charge density, and current 
density. The relative weighting of the contributions is 
controlled by the relative sizes of the components of b^. 
In fact, the same is also true of contributions to gravit- 
omagnetic and higher-order gravitational effects arising 
from other components of the gravitational equations. 
In the interpretation as bumblebee electrodynamics, the 
equations show that we can expect deviations from the 
usual Newton and Coulomb force laws when (3 is nonzero. 
Since the effects are controlled by components of 6 M , they 
are perceived differently from different observer frames, 
much like charge and current in ordinary electrodynam- 
ics. Also, the unconventional source terms are linear in 
f3, so the solutions can be interpreted as superpositions of 
conventional gravitational and electrodynamic fields with 
effects due to (3. However, since (3 itself is formed as the 
combination of the metric and bumblebee fields, the 
extra source terms in fact reflect contributions from the 
gravitational and electrodynamic fields that are absent 
in the Einstein-Maxwell limit. 



3. Timelike case 

In this subsection, we consider in more detail the weak 
static limit with timelike 6 M . We solve explicitly the weak 
static equations in the special observer frame for which 
t>i_i = {b, 0, 0, 0) is purely timelike, and we provide some 
remarks about the more general case. The special frame 
might be identified with the natural frame of the cosmic 
microwave background. However, in laboratory searches 
for modifications to weak static gravitation and electro- 
dynamics, this frame can at best be an approximation 
because the Earth rotates and revolves about the Sun, 
which in turn is moving with respect to the cosmic mi- 
crowave background. Nonetheless, the explicit solutions 
obtained here offer useful insight into the modifications to 
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the Newton and Coulomb potentials that are introduced 
by the massive mode. 

For the purely timelike case, both and B^ v are ex- 
pected to be time independent in the weak static limit. 
However, time dependence must appear in £ even in the 
static limit because as (3 — > the match to electrodynam- 
ics arises in the temporal gauge, for which the standard 
Coulomb solution has a linear time dependence. For this 
reason, in taking the weak static limit in what follows, we 
keep time derivatives acting on the vector components £ . 

The analysis for purely timelike could be performed 
in the harmonic gauge, but it turns out that choosing the 
transverse gauge instead results in the complete decou- 
pling from the massive mode of nine of the 10 gravita- 
tional equations. The transverse gauge is defined by 



d J h i = 0, d J h jk = \d k h 3 -, 



(107) 



loj = u, wrijk — 3-„ . 

where j, k range over the spatial directions. In this gauge 
> g , where 



Goo ~ — V 2 /iqo = 2V 2 $ 9 , where <I> S is the gravitational 



oo 

potential. 

Taking the weak static limit of the gravitational and 
bumblebee equations of motion in this gauge yields 



V 2 $ fi 



V 2 £ - V • d £ 
UE - V(V-f) 



4 7 rG(T°° + T°° 1 -4 K & 3 /?), 
p q - Anb 2 l3, 

J. (108) 



The other nine components of h^ v obey equations that 
are identical to those of general relativity in the weak 
static limit. In Eq. (|108p . is the static energy density 
in the matter sector, and is the static 00 component 
of the zero-/? limit of given in Eq. (|84[) . The charge 
density p q and current density J are also understood to 
be time-independent distributions. 

The explicit form of the massive mode (3 follows from 
Eq. |68|) and is found to be 



p = £o ~ b$ g . 
For purely timelike £> M , the constraint 
d P » 0, 



(109) 



reduces to 



(110) 



so [3 is time independent at this order. As described in 
the previous subsection, (3 can be interpreted as a source 
for $ g and Here, (3 acts as a static source of energy 
density and charge density with 



= bpBV^V 



Op. 0i> 



J% = (p B , 0,0,0), 



(111) 



where pb = —^Kb 2 j3. 

To find the explicit modifications to the gravitational 
and electromagnetic potentials caused by (3, we focus on 
the case of a point particle of mass m and charge q located 
at the origin. The source terms in this case have the form 



and all other components of and are zero. In the 
absence of Lorentz violation, the weak static solution to 
the Einstein-Maxwell equations with these source terms 
consists of the linearized Reissner-Nordstrom metric and 
corresponding electromagnetic fields. Denote by $ m the 
associated gravitational potential. It obeys 



V 2 $ m = 47rG(p„ 



rpOO 



(113) 



In the limit q — > 0, it reduces to the Newton gravita- 
tional potential and obeys the usual Poisson equation 
V 2 <I> m = 4irGp m . Similarly, denote the conventional 
Coulomb potential by $ 9 . It obeys the Poisson equation 



V 2 $ 



(114) 



It is also convenient to express the massive mode (3 in 
terms of a potential $ b ■ We define 







rV 2 $ B (af), 



(115) 



where \Mp\ 2 is the absolute value of the squared mass in 
Eq. (|9"9"|) . Note that <!>b is time- independent due to the 
constraint (|110[) . Since (3 acts as a source of static charge 
density given by Eq. (jllip . the definition of $b can be 
interpreted as a third Poisson equation, 



V 2 $; 



-Pb- 



(116) 



The time independence of both [3 and <&b means they 
are determined by their initial values. 

For the point particle of mass m and charge q, the 
general solution for the gravitational potential <i> g and 
the bumblebee field E^ in the weak static limit can be 
expressed in terms of the conventional potentials $ m and 
$ g and the bumblebee potential $b- We obtain 



* 9 = 
Eo - 



$ m - 4irGb<P B , 
b<f> m -4nGb 2 <I>B 



V 2 $ B , 



Ej = td j [$ q + b<f> r , 



(1 - 47rG6 2 )$ B 



rV 2 $ B ]. 
(117) 



The corresponding static gravitational field g, static elec- 
tric field E, and static magnetic field B are found to be 



g = -V$ m +4ttG6V$ b , 
E 

B = 0. 



-V$ 9 - V$ B , 



(118) 



The static Maxwell equations include a modified Gauss 
law 



rpOO 

1 M 
J M 



Pm 
Pq 



qS^(x), 



(112) 



V-E = -V 2 $ 9 -V 2 $ B 

= Pq+ PB, 



(119) 
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together with the usual static laws 

V-B = 0, 
Vx£ = 0, 
VxB = 0. 



One could, for example, consider initial conditions with 
$b proportional to $ m , 



(120) 



These equations demonstrate for the purely timelike case 
the ways in which a nonzero massive mode (3 modifies 
the conventional static gravitational and electrodynamic 
fields. 

Compared to electrodynamics, the purely timelike KS 
bumblebee model introduces an additional degree of free- 
dom into the problem of determining the static electric 
field E. This extra degree of freedom is the massive mode 
/3, and the extent of its effects depends on the initial con- 
ditions for j3. In the absence of a satisfactory underlying 
theory predicting (3 or of direct experimental observation 
of f3, these initial conditions are undetermined. 

One natural choice is to adopt §b = as an initial 
condition, which implies <&b{x) ~ for all time and 
hence corresponds to zero massive mode (3. The solution 
(|117| then reduces to the weak static limit of Einstein- 
Maxwell theory, as expected. The gravitational poten- 
tial becomes $ ff = <3? m . The bumblebee field reduces to 
£ M = (b& m , tdj(& q +b$> m )), which despite the appearance 
of the gravitational potential implies E = — V$ 9 and the 
usual Gauss and Poisson equations V- E = — V 2< I> g = p q . 
We note in passing that conventional electrodynamics is 
also recovered in the limit l-ftf^ 2 — > oo because taking 
this limit in Eq. (|108[) implies (3 — > and hence — > 0. 
Thus, even with a nonzero massive mode f3, a theory 
with a large mass parameter \Mp \ approximates Einstein- 
Maxwell theory in the low-energy regime, in agreement 
with the discussion following Eq. ([M]) . 

Other choices of initial conditions on $b might appear 
in the context of a more fundamental theory for which a 
bumblebee model provides an effective partial low-energy 
theory. One set of examples with nonzero consists of 
solutions with zero matter current, = 0. Imposing 
Jy, = by hand is a common approach in the litera- 
ture. It implies the bumblebee field has no couplings to 
matter and hence is unrelated to electrodynamics, being 
a field that interacts only gravitationally. The massless 
Lorentz NG modes then propagate freely as sterile fields 
that carry energy and momentum but convey no direct 
forces on charged particles. We emphasize here that a 
theory of this type can nonetheless lead to modifications 
of the gravitational interaction, contrary to some claims 
in the literature. This is readily illustrated in the context 
of the purely timclikc example considered in this subsec- 
tion. With = 0, the linearized solutions for a static 
point particle of mass m are still given by Eqs. (|117[) . 
but Q q — and the fields E and B are largely irrelevant. 
However, the gravitational potential $ 9 is modified by 
the potential <&b for the massive mode /3, with the value 
of $b fixed by initial conditions. 

In the absence of additional theoretical or experimental 
information, the possibilities for model building are vast. 



$i 



a 



4nGb 



(121) 



with a a constant. For this class of solutions, $ g has the 
usual form $ g = —G*m/r for a point mass but with a 
scaled value of the Newton gravitational constant 



G* = (l-a)G N , 



(122) 



where Gat is the value of G* for &b = 0. This rescaling 
has no observable effects in the special observer frame 
because laboratory measurements would determine G* 
instead of Gn, although sidereal, annual, and solar mo- 
tions might introduce detectable effects in realistic exper- 
iments. One could equally well consider instead examples 
in which $s has nonlinear functional dependence on $ m 
or has other coordinate dependence. Various forms could 
be proposed as candidate models to explain phenomena 
such as dark matter or possibly dark energy. With the 
choice of $s being conjecture at the level of the lin- 
earized equations, this approach is purely phenomeno- 
logical. However, if a bumblebee model appears as part 
of the effective theory in a complete theory of quantum 
gravity, the form of $s could well be fixed and predic- 
tions for dark matter and perhaps dark energy could be- 
come possible. 

We conclude this subsection with some comments 
about the case where 6 M takes a more general timelike 
form in the given observer frame O. To investigate this, 
one can either study directly the weak static limit of the 
gravitational and bumblebee equations with a charged 
massive point particle at rest at the origin in O as before, 
or one can boost the observer frame O to another frame 
O' in which 6 M is purely timelike but the point particle is 
moving. The two pictures are related by observer Lorentz 
transformations and are therefore physically equivalent. 

In the frame O in which the particle is at rest, the 
matter energy-momentum tensor and current are given 
by Eq. (|112p as before, but the bumblebee energy- 
momentum tensor and current have additional 
nonzero components due to the nonzero spacelike com- 
ponents of b^. The complete solution therefore involves 
components of /i M „ other than the gravitational poten- 
tial $ g . However, we can gain insight into the behavior 
of the massive mode by considering the constraint (|89")1 
and performing an observer Lorentz boost with velocity 
vb = —b/bo to a frame O' in which 6 M becomes purely 
timelike. Since this transformation must reduce the con- 
straint (|59")l to the purely timelike form (|110|) . it follows 
that in the original frame O the constraint takes the form 



which implies 



d p + v B -V(3 = 0, 



/3 = /3(x-v B t). 



(123) 



(124) 
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We see that the massive mode f3 moves in the frame O 
with a velocity vb equal to the relative boost velocity 
linking to its purely timelike form. The gravitational 
and electrodynamic equations are therefore modified by 
time-dependent source terms moving with velocity vb 
relative to the point particle. An equivalent result is 
obtained by working directly in the frame O' , in which 
is purely timelike and (3 is therefore static, but in which 
all components of the matter energy-momentum tensor 
and current are nonzero and represent a particle moving 
with velocity — vb relative to (3. As before, the form of 
the massive mode (3 is set by initial conditions. Since 
the equations are linear in either frame, the solutions 
consist of a superposition of conventional potentials and 
the massive-mode contributions with relative motion. 



C. Linear Lagrange-multiplier potential 

In this subsection, we discuss the KS bumblebee model 
([26]) with the linear Lagrange-multiplier potential ([6]) and 
the definition (l22l): 



V = V L (X,X) = X(B^B u ±b 2 ). 



(125) 



Bumblebee models with this potential have been widely 
studied at the linearized level for about two decades Q. 
Here, we compare and contrast this theory to the case 
with the smooth quadratic potential Vs- We show the 
Lagrange multiplier A produces effects in the Vl model 
that are very similar to those of the massive mode (3 in 
the Vs model. Paralleling the treatment of the Vs model 
in Sec. IIVB1 the mode expansion = 6 M + £ M of Eq. 
(p0"|) in a Minkowski background is adopted, and the fields 
h^u and £^ are assumed weak so that linearization can 
be performed. 

Variation of the action produces the gravitational and 
bumblebee equations of motion (|72[) and (|77|) . together 
with the Lagrange-multiplier constraint 



X = B^ V B V ± b 2 = 0. 



(126) 



This condition enforces the vanishing of the massive 
mode, as discussed in Sec. Ill Bt which leaves only the 
Lorentz NG modes as possible excitations of the bumble- 
bee field B^. However, there is also an additional degree 
of freedom, the Lagrange multiplier A itself, that appears 
in the equations of motion. 

One might naively expect the Vl model to yield so- 
lutions identical to those obtained in the infinite-mass 
limit \Mp\ — > oo of the Vs model. For example, one 
might reason that an infinite mass would make energeti- 
cally impossible any field excitations away from the mini- 
mum of Vs, leading to the constraint Q126p . However, the 
potential Vl is a function of two combinations of fields, 
Vl = Vl(X,X), whereas Vs(X) involves only one. The 
infinite-mass limit indeed suppresses X excitations away 
from X = in Vs, but it contains no match for A. There 
is thus an extra field degree of freedom in Vl relative to 



the infinite-mass limit of Vs- For example, Vg — > in 
the infinite-mass limit because no excitations of X are 
allowed, whereas V[ = A, which need not vanish. We 
therefore conclude that the correspondence between the 
infinite-mass limit of the Vs model and the Vl model can 
occur only when A = 0. 

To gain intuition about the effects associated with the 
Lagrange-multiplier field A, consider its role in the equa- 
tions of motion. Since the constraint (|126|) ensures no 
massive excitations can appear, the bumblebee energy- 
momentum tensor (|74l) reduces to 



1 B 



EM 



J y ■ 



We find 



1 V 



2XB^B U « 2¥-b v X, 



(127) 



(128) 



where the last form is the leading-order contribution in 
the linearized limit. Similarly, the bumblebee current 
(|75|) becomes 



(129) 



Conservation of the matter current J*,, implies the con- 



straint 



b^d^X « 0. 



(130) 



When A — > 0, all these equations reduce to those of 
Einstein-Maxwell theory, in agreement with the discus- 
sion of bumblebee electrodynamics in Sec. IIV A 31 More- 
over, by comparison with Eqs. ([57)) - (|89"|) , we see that the 
Lagrange multiplier A in the Vl model plays a role very 
comparable to that of the massive mode [3 in the Vs 
model. In effect, A acts as an additional source of energy- 
momentum density and current density in the equations 
of motion. 

The propagating modes for the KS bumblebee model 
with the linear Lagrange-multiplier potential Vl have 
been investigated elsewhere [!, Q • The usual two gravi- 
ton modes propagate as free transverse massless modes 
independently of A, as do the usual two photon modes 
emerging from the Lorentz NG modes. Since A has no 
kinetic terms, it is auxiliary and cannot propagate. 

The weak static limit of the Vl model can also be stud- 
ied, following the lines of the discussion for the Vs model 
in Sec. IIV B 31 Consider, for example, the case of purely 
timelike 6 M = (6,0,0,0). The constraint (|130[) then im- 
plies that A must be time independent at leading order, 



dnX « 0. 



(131) 



A nonzero A therefore acts as an additional static source 
of energy density and charge density, which can mod- 
ify the static potentials. We again adopt the trans- 
verse gauge (|107p . for which the gravitational potential 
is $ g = — hoo/2. Linearizing the constraint (|126[) for the 
purely timclike case yields 



£ - = 0. 



(132) 
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As expected, this condition corresponds in the context of 
the Vs model to enforcing the vanishing of the massive 
mode p. In the context of Einstein-Maxwell theory, it 
corresponds as before to a gauge-fixing condition that 
reduces to the usual temporal gauge in the absence of 
gravity. 

For a single point particle of charge q and mass m at 
rest at the origin, the energy density p m and charge den- 
sity p q are given by Eq. (|112p . The equations of motion 
in the weak static limit become 

V 2 $ 3 « 4 7 rG(p m +T™ 1 + 2& 2 A), 
V 2 £ -V-d 8 w p q + 2bX, 
□ f-V(V-f) « 0. (133) 

Comparison of these expressions with Eqs. (|108p for the 
Vs model again reveals the correspondence between the 
roles of the Lagrange-multiplier field A and the massive 
mode (5. 

With A acting as an effective source of charge and en- 
ergy density, we can introduce a potential <f>s, 

A = -^V 2 $ B , (134) 

which cquivalently can be viewed as the solution to the 
Poisson equation 

V 2 $ B = -p B , p B = 2b\. (135) 

As before, let $ m denote the usual gravitational potential 
obeying Eq. (] 1 1 3[) . and let <j> g denote the usual Coulomb 
potential obeying Eq. (|114f> . In terms of these three po- 
tentials, the weak-field static solutions to the equations 
of motion (| 1 33[) can be written as 

$ g = $ m -47rG&$B, 
So = 6$ m - 4vrGfe 2 $B, 

£ 3 = td 3 [<f> q + b<S> m + {l-4irGb 2 )<5> B }- (136) 

The associated static gravitational field g, static electric 
field E, and static magnetic field B are all given by the 
same mathematical expressions as Eqs. (Ill 8[) . but the 
potential $b is now defined in terms of A according to 
Eq. (|134[) . Similarly, a nonzero Lagrange- multipler field 
A yields the same modified mathematical form (|119[) of 
the Gauss law, but the potential and charge pb are 
now defined in terms of A. 

The bumblebee potential $b in Eqs. (|136p is unde- 
termined by the equations of motion (|133p and must be 
specified as an initial value. This specification also fixes 
the initial value of the Lagrange multiplier A, and the 
condition (|131[) then ensures that A remains unchanged 
for all time. The situation here parallels that of the Vs 
model, where initial conditions must be imposed that 
subsequently fix the bumblebee potential and the massive 
mode (3 for all time. These results are special cases of a 
more general fact often overlooked in the literature: to be 
well defined, all bumblebee models require explicit initial 



conditions on the massive modes and Lagrange-multiplier 
fields. The subsequent development of the massive modes 
and Lagrange multipliers can then be deduced from the 
equations of motion or from derived constraints. In the 
absence of specified initial conditions, physical interpre- 
tations and predictions cannot be reliable. Moreover, in 
the absence of direct experimental observation or pre- 
diction from an underlying theory, the choice of initial 
conditions is largely unrestricted and can lead to widely 
differing effects. 

One natural initial condition is $b = 0, A = 0. The 
above results then reduce to Einstein-Maxwell theory in 
the weak static limit, as expected. This Vl model has 
only field excitations maintaining both Vl — and V' h — 
0. It corresponds to the infinite-mass limit of Vs, which 
is itself equivalent to Einstein-Maxwell theory. 

Another possibility is to consider models without di- 
rect coupling to matter, Jj^ = 0, which implies p q = 0, 
$5 = in the above equations. Like their Vs counter- 
parts discussed in Sec. IIV B 3[ these models are purely 
gravitational and are unrelated to electrodynamics. As 
always, initial conditions on <5>b must be specified. One 
simple possibility, for example, is to choose $b = — bjltq- 
This Vl model is contained in the analysis of Ref. [471 ] . 
The spatial components of the bumblebee vanish, £j = 0, 
so the bumblebee field B^ is parallel to a timelike Killing 
vector. The modifications to $ g in this example involve 
a rescaling of the Newton gravitational constant. 

Examples with both and $ 9 nonzero can also be 
considered. These incorporate direct charge couplings 
to matter, but the static weak-field limit yields modi- 
fied gravitational and Coulomb potentials given by Eq. 
(|136[) . One simple example is the choice $s = — $ g . 
Since <E> g ^ $ m , this Vl model has a modified gravita- 
tional potential. The solution for the bumblebee field 
has the form of a total derivative, £^ = d^itb^g). In 
the limit q — > and A — > 0, it provides an example of 
a solution that with A = and hence $b — has been 
identified in Ref. [45| as potentially flawed due to the for- 
mation of shock discontinuities in £^. Whether or not A 
is zero, this solution is unusual. Both the field strength 
and the energy-momentum tensor vanish be- 
cause the effective charge density ps associated with the 
Lagrange-multiplier field cancels the matter charge den- 
sity p q . Although the field strengths are zero, the bum- 
blebee field £ il must be nonzero because the constraint 
(|132[) implies £q — b$ g , which cannot vanish in the pres- 
ence of gravity. For a point charge the solution £ fJi does 
indeed contain a singularity, but this is physically unre- 
markable as it merely reflects the usual 1/r dependence 
in The same behavior arises in the standard solu- 
tions of Einstein-Maxwell theory in a gauge fixed by Eq. 
(|132p . to which the Vl model is equivalent. 
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D. Quadratic Lagrange-multiplier potential 

As a final example, we consider in this subsection the 
specific KS bumblebee model (I26p with the quadratic 
Lagrange-multiplier potential ([7]) 

V = V Q {\X)= 1 1 \(B,g^B v ±b 2 ) 2 , (137) 

where X is defined in Eq. (|22|) . As before, we adopt the 
mode expansion = b^+E^ of Eq. (|60p in a Minkowski 
background, and where useful we assume weak fields /i M „ 
and £ M . 

Bumblebee models with a quadratic Lagrange- 
multiplier potential have not previously been considered 
in detail in the literature. We introduce them here partly 
as a foil for the Vl case, in which the Lagrange mul- 
tiplier plays a key role in the physics of the model. In 
contrast, the Lagrange multiplier for the potential Vq de- 
couples entirely from the classical dynamics. The point 
is that variation with respect to A yields the constraint 
X 2 = 0, which is equivalent to X — and forces the mas- 
sive mode to vanish. However, the quadratic nature of 
Vq means that A always appears multiplied by X in the 
gravitational and bumblebee field equations, so the on- 
shell condition X = forces the field A to decouple. Also, 
the potential obeys both Vq = and Vq — on shell, so 
the bumblebee energy-momentum tensor and bumblebee 
current reduce to the standard expressions for electrody- 
namics in curved spacetime. The equations of motion are 
therefore equivalent to those of Einstein-Maxwell electro- 
dynamics in the gauge X = 0. This correspondence holds 
both in the linearized limit and in the full nonlinear the- 
ory. 

Since the equations of motion generated by a theory 
with a quadratic Lagrange-multiplier potential are the 
equations in the absence of the constraint plus the con- 
straint itself, introducing this type of potential is equiv- 
alent at the classical level to imposing the constraint by 
hand on the equations of motion. The only distinction 
between the two approaches is the presence of the decou- 
pled Lagrange-multiplier field. Bumblebee models with 
the potential Vq therefore incorporate a variety of models 
in which constraints are imposed by hand. For example, 
Will and Nordtvedt have considered solutions involving 
a nonzero background value for a vector in bumblebee- 
type models with Lagrange density of the form (fTi?)) but 
without a potential term [42l |. This approach produces 
equations of motion identical to those of the correspond- 
ing Vq model. 

The Vq models are also related to theories in which the 
constraint is substituted into the action before varying, 
although the correspondence is inexact. For example, 
Nambu has investigated the Maxwell action in Minkowski 
spacetime using the constraint X = for the purely time- 
like case as a nonlinear gauge condition substituted into 
the action prior to the variation [53|. This represents 
gauge fixing at the level of the action, and it yields a to- 
tal of four equations for the fields: the original constraint 



and three equations of motion from the variation. How- 
ever, the corresponding Vq model yields five equations 
of motion instead, one of which is the constraint. The 
extra equation is the Gauss law, which in the Nambu 
approach is imposed as a separate initial condition that 
subsequently holds at all times by virtue of the three 
equations of motion and the constraint. 

For the specific KS bumblebee model (f2l)]) with poten- 
tial Vq in Eq. (|137[) . the freely propagating modes can 
be found by linearizing the gravitational and bumble- 
bee equations of motion (f?2")) and (|77|) and the constraint 
X = 0. The linearization generates the same equations 
as emerge for the Vl model with potential (|125[) in the 
limit A — > 0. The propagating modes consist of the usual 
graviton and photon in an axial gauge as expected. 
Similarly, the weak static limit of the Vq model with po- 
tential (|137|) produces equations identical at linear level 
to those of the Vl model with potential (|125[) in the 
limit A — > 0. The usual gravitational and electromagnetic 
Poisson equations therefore emerge, and the correct New- 
ton and Coulomb potentials hold at the linearized level. 

The exact correspondence between the static limit of 
the Vq bumblebee model and the Newton and Coulomb 
potentials involves the nonlinear constraint (| 1 26[) . The 
explicit forms of the solutions are therefore also nonlin- 
ear. For example, consider again the case of a point 
particle of mass m and charge q at rest at the ori- 
gin in the presence of a purely timelike vacuum value 
= (6,0,0,0). The equation of motion for A is the 
constraint (|126p . which represents a nonlinear condition 
relating £^ and h^ v . Expanding this constraint through 
quadratic order gives 

£ 2 + 2b£ - 26 2 $ g - ib$ m £ - £ 2 « 0, (138) 

where $ m obeys the conventional Poisson equation Q113P 
for the point source with mass and charge density given 
by Eq. pi2[) . In the transverse gauge, the gravitational 
potential Q g obeys 

V 2 $ ff «47rG(p m +T^), (139) 

where the energy-momentum tensor has the usual 
Einstein-Maxwell form . The solution for is there- 
fore the standard gravitational potential for a static point 
charge in a curved spacetime in the weak-field limit. At 
quadratic order, the solution £^ — (£o,£j) satisfying the 
constraint (|138[) and the field equations is found to be 

36 t 2 

£ « b<P g + -<£ 2 g + -{d J (<P q + b<P g )} 2 , 
£j « td j ($ q + b$ g ) + 3bt$ g d j $ g 

+^[9 fc ($ g + bQ^djdkiQg + b$ g ), (140) 

where $ 9 is the conventional Coulomb potential obeying 
the Poisson equation (| 1 14[) . These expressions are both 
time dependent and nonlinear, but they nonetheless gen- 
erate the usual static electric field E and magnetic field 
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B for a point charge. Explicit calculation to quadratic 
order yields 

S«-V$ g , BwO, (141) 

which implies the usual form V • E « p 9 of the Gauss 
law. From these equations, we see again that Einstein- 
Maxwell theory is recovered despite the absence of U(l) 
symmetry. In effect, the nonlinear condition Q138P plays 
the role of a nonlinear gauge-fixing condition in a U(l)- 
invariant theory, removing a degree of freedom and leav- 
ing only NG modes that propagate as photons and gener- 
ate the usual Coulomb potential in the weak-field limit. 



V. SUMMARY 

In this paper, we have investigated the properties of 
the massive modes that can emerge from spontaneous 
local Lorentz and diffcomorphism breaking. In Ricmann 
spacetime, no massive modes of the conventional Higgs 
type can appear because covariant kinetic terms involve 
connections with derivatives. However, an alternative 
form of the Higgs mechanism can occur instead, in which 
massive modes originate from quadratic couplings in the 
potential V inducing the symmetry breaking. 

SectionlTTlprovides an analysis of this alternative Higgs 
mechanism in the general context of an arbitrary tensor 
field. Both smooth potentials and Lagrange-multiplier 
potentials are considered. Massive modes appear for a 
smooth potential when excitations with V' =/= exist, 
and they are formed as combinations of field and metric 
fluctuations. For Lagrange-multiplier potentials the mas- 
sive modes are constrained to vanish, but the Lagrange 
multiplier fields can play a related role. The propagation 
of the massive modes depends on the nature of the ki- 
netic terms in the theory, and the requirements of unitar- 
ity and ghost-free propagation constrain possible models. 
Even if the massive modes are nonpropagating they can 
influence gravitational phenomena through, for example, 
effects on the static gravitational potential. These mod- 
ifications are of potential interest in alternative theories 
of gravity and descriptions of phenomena such as dark 
matter or dark energy. 

Following the general treatment, we investigate in Sec. 
Illll a broad class of theories called bumblebee models that 
involve gravitationally coupled vector fields with sponta- 
neous local Lorentz and diffeomorphism breaking. For 
arbitrary quadratic kinetic terms, the Lagrange density 
is given in Eq. (|19p. Along with the symmetry-breaking 
potential V and a matter sector, this Lagrange density 
involves five parameters, four of which can be linearly 
independent in specific models. A particularly attractive 
class of theories are the KS bumblebee models, which 
have kinetic term of the Maxwell form and hence an addi- 
tional constraint that minimizes problems with unitarity 
and ghosts. These models also offer candidate alterna- 
tives to Einstein-Maxwell electrodynamics. 



In a series of subsections and the associated appendix, 
we provide some results valid for general bumblebee mod- 
els. The observer and particle forms of local Lorentz 
transformations and diffcomorphisms are presented. Us- 
ing the vierbein formalism, some decompositions of the 
bumblebee field and metric are given that are suitable 
for the identification of Lorentz and diffeomorphism NG 
modes. The effects of various choices of Lorentz and dif- 
feomorphism gauges are described. Alternative decom- 
positions used in some of the literature are also discussed, 
in which the Lorentz NG modes are hidden and only 
spacetime variables are used. 

To provide explicit examples and to gain insight via a 
more detailed analysis, we focus attention in Sec. IIVI on 
the class of KS bumblebee models. The basic equations 
of motion and conservation laws are obtained, and some 
properties of the bumblebee currents are considered. The 
interpretation of these bumblebee models as theories of 
electromagnetism and gravity, known as bumblebee elec- 
trodynamics, is discussed. They contain four transverse 
massless modes, two of which are massless gravitons and 
two of which are massless photons, along with a massive 
mode or Lagrange-multiplier mode. When the massive 
mode or Lagrange multiplcr vanishes, or in the limit of 
infinite mass, conventional Einstein-Maxwell theory in an 
axial gauge is recovered. 

Section IIVI also contains subsections considering in 
more detail various types of potentials, including the 
smooth potential Vs in Eq. the linear Lagrange- 

multiplier potential Vl in Eq. (|125[) . and the quadratic 
Lagrange-multiplier potential Vq in Eq. (|137p . For the 
Vs model, the gravitational and bumblebee equations of 
motion are investigated to determine whether a physical 
massive mode can propagate as a free field. The sign of 
the squared mass term depends on whether the bumble- 
bee vacuum value is timelike or spacelike. In the time- 
like case, the massive mode is a ghost, while in the space- 
like case the squared mass has the usual sign. However, 
in both cases the dispersion law is unconventional and 
no localized physical solutions satisfying suitable asymp- 
totic boundary conditions exist, so the Vs model has no 
freely propagating massive modes. Nonetheless, the mas- 
sive mode has a physical impact as an auxiliary field, 
acting as an additional source of energy-momentum den- 
sity and current density. In the weak static limit, for 
example, the solutions to the equations of motion in the 
presence of mass and charge describe modified Newton 
and Coulomb potentials according to Eq. (|117p . These 
may be of phenomcnological interest in the context of 
dark matter and perhaps also dark energy. The effects 
are controlled by the massive mode, but its form is dy- 
namically undetermined and must be imposed via initial 
conditions. 

Many of the results obtained for the Vs model ap- 
ply also to the Vl and Vq models that have Lagrange- 
multiplier fields. For example, all the models contain four 
massless propagating modes behaving like gravitons and 
photons. Although the Vl and Vq models generate an 
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additional constraint that eliminates the massive mode 
on shell, the Lagrange-multiplier field appears instead as 
a extra degree of freedom playing a similar role to that of 
the massive mode in the Vs model. The form of the La- 
grange multiplier must be set by initial conditions, and 
if it is chosen to vanish then Einstein-Maxwell theory is 
recovered. The key difference between the Vl and Vq 
models lies in the role of the Lagrange multiplier, which 
can affect the physics as an auxiliary mode in the Vl 
model but which decouples from the theory in the Vq 
model. 

In the context of bumblebee electrodynamics, the mas- 
sive mode or Lagrange-multiplier field acts as an ad- 
ditional degree of freedom relative to Einstein-Maxwell 
theory. The extra freedom arises because the bumblebee 
model has no U(l) gauge symmetry, and the structure 
of the kinetic terms implies that the freedom must be 
fixed as an initial condition. In the absence of experi- 
mental evidence for a massive mode or of guidance from 
an underlying theory, the choice of initial condition is 
largely arbitrary. A natural choice sets the massive mode 
or Lagrange multiplier field to zero, reducing the the- 
ory to Einstein-Maxwell electrodynamics up to possible 
SME matter-sector couplings. Bumblebee electrodynam- 
ics therefore provides a candidate alternative explanation 
for the existence of massless photons, based on the mass- 
lessness of NG modes instead of the usual gauge symme- 
try. In any case, the possibility of spontaneous breaking 
of local Lorentz and diffeomorphism symmetry remains 
a promising avenue for exploring physics emerging from 
the Planck scale. 



Acknowledgments 

This work was supported in part by DOE grant DE- 
FG02-91ER40661, NASA grant NAG3-2914, and NSF 
grant PHY-0554663. 



APPENDIX A: TRANSFORMATIONS 

In this appendix, we provide explicit transformation 
formulae for the vierbein, bumblebee, and NG field 
components under particle diffeomorphisms, particle lo- 
cal Lorentz transformations, observer general coordinate 
transformations, and observer local Lorentz transforma- 
tions. The vierbein decomposition is given in Eq. (|37[) . 
the bumblebee decomposition is given in Eq. (|40[) . and 
the NG fields are defined in Eq. (j4"T|) . The transformation 
formulae can be deduced from these expressions and from 
the behavior of the vacuum expectation values. All the 
formulae given here are linearized assuming infinitesimal 
fields and transformation parameters. 



1. Particle transformations 

Under infinitesimal particle diffeomorphisms, 
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Under infinitesimal local particle Lorentz transforma- 
tions, 
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2. Observer transformations 

The primary difference between observer and parti- 
cle transformations is that the components of vacuum- 
valued fields transform under observer transformations. 
However, additional differences arise when the distinc- 
tion between spacetime and local indices is dropped. In 
particular, care is required in expressing observer trans- 
formation laws for the metric and the vierbein. 

Three versions of the Minkowski metric appear in the 
formalism: the spacetime metric rj^ u , the local metric 
r) a b, and the mixed metric rj^a used with the vierbein. All 
three are numerically equal in cartesian coordinates, but 
they behave differently under observer transformations. 
In the customary notation using only Greek indices, these 
three metrics must be labeled differently. Since they rep- 
resent vacuum values of the metric, we write them as 
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expectation values and define 



eral coordinate transformations are 
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The inverse metrics are denoted with upper indices and 
are defined so that 



(A4) 



In terms of these vacuum values, the metric and its in- 
verse involve infinitesimal excitations: 



(A5) 



The transformation laws for the vierbein also depend 
on the component basis. We define 
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For example, these vacuum values obey 
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as expected from the relationship between the metric and 
the vierbein. The vierbein field itself can be written 



&p,v ~ (j)iiv) e "i" (Xftf) "i" 2 Xftf! 

e""« (^)e + (x^)-^^+X AW . 



(A8) 



In much of the literature, the vacuum value for (x^v) 
is assumed to vanish in cartesian coordinates. How- 
ever, this choice is observer dependent. Under cither 
observer general coordinate transformations or observer 
local Lorentz transformations, a zero value of (x M „) trans- 
forms into a nonzero value. We therefore include the 
vacuum value (Xfw) in the formulae below. 

The transformation rules for infinitesimal observer gen- 
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The formulae for infinitesimal observer local Lorentz 
transformations are 
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In these expressions, indices can be raised and lowered 
using (f]^ u ) n . 

For some purposes, it may be desirable for an ob- 
server transformation to take the appearance of the in- 
verse of a particle transformation. This can be achieved 
by changing the signs of the parameters, £^ — ► — ^ and 
e uv — > — e uv , cither in the formulae for the observer 
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transformations or in those for the particle transforma- 
tions. For example, with signs chosen appropriately the 
bumblebee field i? M can be arranged to be invariant, 
— > , under the composite transformation consisting 
of a particle diffeomorphism followed by the correspond- 
ing inverse general coordinate transformation. However, 



the reader is cautioned that invariance of this type may 
fail for the vacuum value and excitations of B^. For in- 
stance, the vacuum value is invariant under particle 
diffcomorphisms but transforms as ordinary components 
of a vector under observer general coordinate transfor- 
mations. 
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